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ABSTRACT 

'The  lack  of  reliable  data  on  the  fluidization  and  heat  transfer 
characteristics  of  low  density  particles  in  a  fluidized  bed  has 
prompted  an  experimental  and  analytical  investigation  into  this 
subject.  Seven  groups  of  particles  ranging  in  diameter  from  0.25  mm 
to  2.0  mm  and  density  from  2.5  to  32  pcf  have  been  successfully 
fluidized  and  have  t>eeri; shown  to  be  generally  well  predicted  by 
classical  fluidization  and  fluidized  bed  heat  transfer  theory.  Two 
other  groups  of  particles,  also  in  this  approximate  range  of  particle 
diameter  and  density,  are,  however,  unable  to  be  fluidized  due  to 
significant  inter-particle  and  static  electric  attractions. 

Using  the  experimental  data  and  results  as  a  basis  of  analysis, 
two  application  of  low  density  particle  fluidization  in  a  building 
efficient  energy  management  program  are  discussed.  A  fluidized  bed 
can  be  incorporated  into  the  wall  cavity  of  a  building  for  use  as 
either  a  collector  of  solar  energy  or  as  a  heat  exchange  medium  in  a 
building  space  heating/cooling  program.  As  a  solar  collector,  it  is 
shown  that  the  low  density  particle  fluidized  bed  would  thermally 
perform  between  comparable  conventional  liquid  and  air-cooled  flat 
plate  solar  collectors.  It  would  require  less  water  pumping  power  and 
plumbing  than  the  liquid  collector  and  less  air  pumping  power  than  the 
air  collector.  -  -  •  .  ■  r>  ' 

As  a  heat  exchanger,  the  low  density  particle  fluidized  bed  would 
function  in  a  building  space  comfort  control  capacity  by  providing  a 
large  wall  surface  area  available  for  heat  transfer  to  or  from  the 
space.  It  is  shown  that  this  rate  of  heat  transfer  provided  by  the 
fluidized  bed  is  significantly  greater  than  what  is  available  for  air 
alone  flowing  over  the  surface  and  at  reduced  air  pumping  power 
requirements. 
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NOMENCLATURE 

solar  collector  surface  area 
total  fin  surface  area 

total  inside  surface  area  of  a  heat  transfer  tube 

total  unfinned  surface  area 

cross-sectional  area 

particle  diameter 

gravitational  constant 

fin  heat  transfer  coefficient 

working  fluid  heat  transfer  coefficient  inside  a  tube 
unfinned  heat  transfer  coefficients 
solar  insolation  on  a  collector  surface 
denotes  minimum  fluidization  conditions 

fluidized  bed  pressure  drop  per  unit  height  of  bed 
Prandtl  number  of  the  fluidizing  gas 
rate  of  heat  transfer 

rate  of  heat  loss  from  a  solar  collector  plate 
to  the  environment 

rate  of  heat  transfer  from  the  solar  collector 
plate  to  the  working  fluid  (useful  energy) 

the  temperature  difference  between  a  heated 
surface  and  the  fluidized  bed 


< 


U 


fluidized  bed  bulk  temperature 
heated  surface  temperature 
fluidizing  gas  velocity  at  distributor 


solar  absorbance  of  a  collector  plate  surface 

fraction  of  voids  in  a  bed  of  particles 

emmissivity  of  a  heated  surface 

fin  efficiency 

thermal  conductivity  of  air 

thermal  conductivity  of  a  fluidized  particle 

dynamic  viscosity  of  the  fluidizing  gas 

kinematic  viscosity  of  the  fluidizing  gas 

density  of  the  fluidizing  gas 

density  of  the  fluidized  particles 

Stephan-Boltzman  constant 

effective  solar  transmittance  of  the  solar 
collector  cover(s) 

particle  sphericity 


CHAPTER  1 


INTRODUCTION 


Fluidized  bed  technology  has  been  successfully  applied  to  various 
industrial  process  problems  such  as  catalytic  cracking  and  the  removal  of 
sulfur  based  pollutants  associated  with  coal  combustion.  These  successes 
have  continued  to  spur  research  efforts  into  many  other  areas.  One  area, 
though,  that  has  not  been  explored  until  now  is  the  application  of  air 
fluidized  bed  technology  in  building  energy  management. 

A  review  of  the  literature  indicates  that  previous  investigations  of 
fluidized  bed  technology  have  used  as  a  fluidizing  material  particles 
ranging  in  size  from  approximately  50  um  to  1  mm  and  densities  from 

3 

500  -  2500  kg/m  .  These  materials,  which  generally  include  sand,  metallic 
grit,  glass  beads  and  limestone,  were  selected  on  the  basis  of  being  best 
suited  for  the  intended  applications  and  studies.  Thus,  the  theories, 
practices,  successes  and  failures  of  past  work  in  fluidized  bed  technology 
have  all  been  based  on  the  use  of  these  dense  type  materials. 

As  a  result,  this  extensive  amount  of  work  has  developed  well  known 
relations  and  theories  relating  such  concepts  as  minimum  fluidization 
velocities,  bed  pressure  drops,  and  bed  heat  transfer  to  particle  dimater 
and  density. 

A  goal  of  this  work,  on  the  other  hand,  is  to  assess  the  feasibility 
of  using  a  fluidized  bed  as  part  of  a  building  efficient  energy  management 
program.  The  concept  of  using  a  fluidized  bed  in  this  manner  consists  of 
It  being  integrated  into  a  building  wall  cavity  to  serve  as  a  collector  of 
solar  insulation  or  as  a  heat  exchange  medium. 


In  these  applications  a  building's  wall  cavity  would  be  filled  with  a 
fluidizing  material.  As  a  solar  collector,  the  exterior  side  of  a  south 
facing  wall  would  serve  as  a  solar  insulation  collector  plate.  The  heat 
would  be  transferred  by  the  fluidized  particles  to  immersed  heat  transfer 
tubes.  A  working  fluid,  such  as  water,  passing  through  the  tube  would 
collect  this  energy  and  transport  it  to  another  location  out  of  the  wall 
for  storage  or  immediate  use  (see  Figure  7.1). 

As  a  heat  exchange  medium  the  fluidized  bed  could  be  used  in  either  a 
building  space  heating  or  cooling  scheme.  Due  to  the  fluidized  bed 
isothermal  characteristic,  a  constant  temperature  wall,  either  above  or 
below  room  temperature,  could  be  maintained  to  provide  the  desired  path  of 
heat  transfer. 

For  either  application,  this  use  of  a  fluidized  bed  must  provide 
thermal  performance  and  operational  requirements  that  are  comparable,  or 
better,  to  conventional  alternatives  in  order  to  be  feasible.  Two 
operational  requirements,  bed  pressure  drop  and  external  power  needs, 
mandate  the  use  of  low  density  fluidization  materials  in  these 
applications.  Fluidizing  materials  such  as  sand  or  a  metallic  grit  yield 
large  bed  pressure  drops,  in  excess  of  the  structural  capability  of  a 
typical  building  wall  cavity.  These  materials  also  require  significant  air 
pumping  power  for  fluidization.  Low  density  materials  (less  than  500 
Kg/m  ),  on  the  other  hand,  would  greatly  reduce  these  effects  and  thus 
provide  a  potentially  more  reliable  and  less  expensive  alternative  to 
conventional  equipment. 

However,  a  review  of  the  literature  also  reveals  that  little  data 
exists  on  the  fluidization  and  heat  transfer  of  low  density  materials  in  a 
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fluidized  bed  [1,3,4].  This  evaluation  must  be  conducted  in  order  to  first 
completely  determine  the  feasibility  of  the  use  of  fluidized  beds  of  low 
density  materials  in  a  building  energy  management  program  and  second  to 
provide  a  further  broadening  of  understanding  in  fluidized  bed  technology. 

To  this  end,  the  primary  objective  of  this  work  has  been  to  identify 
and  test  various  low  density  materials  to  determine  their  fluidization  and 
heat  transfer  characteristics.  Fluidization  studies  are  made  to  determine 
minimim  fluidization  velocities,  bed  pressure  drops,  flow  regimes,  and  the 
parameters  that  effect  these.  Heat  transfer  studies  are  conducted  from  a 
heated  bed  wall  surface  {to  simulate  a  heated  solar  collector  plate)  and  an 
immersed  tube  (to  simulate  heating  the  working  fluid).  Data  analysis  is 
based  upon  the  theories  and  trends  previously  developed  with  denser 
materials  to  determine  the  degree  of  correlation  or  discrepancy  between 
experimental  results  and  theoretical  predictions  for  the  low  density 
material. 

Based  upon  the  experimental  results,  comparisons  are  made  between  the 
expressed  applications  of  the  low  density  material  fluidized  bed  and  the 
conventional  counterparts. 

Initially,  however,  a  review  of  the  fundamental  theories  of  fluidized 
bed  technology  must  be  conducted  in  order  to  provide  the  framework  upon 
which  this  investigation  is  based. 


CHAPTER  2 


THEORETICAL  BASIS  OF  INVESTIGATION 


This  chapter  will  outline  the  theory  upon  which  the  experimental 
design  and  analysis  of  this  investigation  is  based.  A  review  of  the 
fundamentals  of  fluidization  and  fluidized  bed  heat  transfer  will  be 
conducted  as  developed  from  previous  work  with  dense  particles.  The  brief 
examination  of  theory  presented  here  will  yield,  in  general  terms,  the 
tools  and  concepts  required  for  the  detailed  investigation  of  low  density 
particles.  The  specific  development  of  these  tools  and  concepts  for  this 
particular  experiment  will  be  conducted  in  Chapter  3. 

2.1  The  Mechanisms  of  Fluidization 

Simply  stated,  fluidization  is  the  operation  by  which  an  agglomeration 
of  solid  particles  are  transformed  into  a  fluid-like  state  when  brought 
into  contact  with  a  fluid.  Although  fluidization  may  be  accomplished  with 
either  a  liquid  or  gas,  the  concern  of  this  investigation  is  fluidization 
with  air  at  ambient  conditions. 

As  the  air  flows  upward  through  a  bed  of  solid  particles,  at  low  flow 
rates,  it  percolates  through  the  void  spaces  between  the  particles.  This 
condition  is  known  as  a  fixed  bed.  However,  as  the  air  flow  rate  increases 
the  particles  begin  to  move  apart  and  may  vibrate  slightly  C73  .  With 
further  increases  of  flow  rate  a  point  is  eventually  reached  where  the  drag 
force  exerted  on  the  individual  particles  by  the  upward  moving  gas  just 
equals  the  weight  of  the  particles,  the  bed  height  has  Increased,  and  the 
particles  have  become  suspended.  This  point  Is  called  minimum 
fluidization.  The  air  velocity  required  to  attain  minimum  fluidization  is 
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called  the  minimum  fluidization  velocity,  ,  and  the  pressure  drop 

within  the  bed  at  this  point  equals  approximatedly  the  cumulative  weight  of 
the  suspended  particles. 

Ergun  [2]  has  developed  a  relationship  for  the  pressure  drop  of  a  gas 
through  a  fixed  bed  of  uniformly  sized  solids  as 
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At  minimum  fluidization,  where  the  pressure  drop  equals  approximately  the 
weight  of  the  suspended  particles  the  following  relation  holds 
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By  combining  Eqs.  (2.1)  and  (2.2)  Kunii  and  Levenspiel  [3]  have 

developed  an  expression  for  analytically  predicting  the  minimum 

fluidization  velocity  of  a  bed  of  uniformly  sized  particles 
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At  air  velocities  beyond  the  minimum  fluidization  velocity,  Uq  >  Um^, 
instabilities  begin  occurring  within  the  bed  in  the  form  of  bubble  creation 
at  the  air  distributor  and  enlargement  as  they  rise  to  the  top  of  the  bed. 
Davidson  [4]  has  attempted  to  explain  this  phenomenon  through  the  use  of 
his  "Two-Phase  Theory  of  Fluidization".  He  states  that  when  U0  >  Umf  a 
fluidized  bed  can  be  visualized  as  having  two  simultaneous  phases:  an 
emulsion  phase  and  a  bubble  phase.  As  air  flow  at  Uq  enters  the  bed,  the 
emulsion  phase  (the  fluidized  particles)  constantly  remains  at  e  f 

with  a  corresponding  Um^  requirement.  The  excess  volume  flow  of  air, 
(Uq  -  Umf)  x  (bed  cross-sectional  area),  passes  through  the  bed  as 
bubbles.  The  presence  of  bubbles  determines  the  total  bed  expansion  and 
provides  for  solid  circulation  within  the  bed.  This  bubbling  condition  is 
known  as  an  aggregatively  fluidized  bed. 

Despite  this  increase  in  air  flow  beyond  minimum  fluidization  and  the 
resultant  bubble  formation  the  bed  pressure  drop  remains  relatively 
unchanged  from  that  at  minimum  fluidization.  This  is  due  to  the  fact  that 
the  dense  gas-solid  phase  is  well  aerated  and  can  deform  easily  without 
appreciable  resistance.  Thus,  the  minimum  fluidization  velocity  can  be 
experimentally  determined  by  varying  Uq  ,  beyond  minimum  fluidization, 
measuring  the  corresponding  bed  pressure  drop,  and  checking  for  its 
constancy  as  Uq  is  decreased  to  the  point  of  minimum  fluidization.  Kunii 
and  Levenspiel  recommend  using  this  experimental  method  of  determining 
U^  for  the  particular  size  distribution  of  particles  actually  in  a  bed 
rather  than  Eq.  (2.3). 

Additionally,  Geldart  [5]  has  classified  materials  into  four  general 
groups  according  to  their  fluidization  characteristics  (Fig.  2.1): 


# 


A.  particles  of  low  density  (less  than  1400  kg/nr)  and 
of  mean  diameter  between  20  and  100  pm; 

B.  powders  having  a  mean  diameter  between  40  and  500 

3 

m  and  a  density  in  the  range  1400  to  4000  kg/m  ; 

C.  small  particles  having  a  mean  size  less  than  30  urn 
for  which  the  effect  of  interparticle  forces  is  significant; 

0.  large  (>600  urn)  and/or  dense  particles. 

Group  A  particles  exhibit  considerable  bed  expansion  before  the  onset  of 
bubbling  within  the  bed  and  collapse  slowly  after  the  air  supply  is  turned 
off.  Materials  in  Group  B  exhibit  much  less  stable  bed  expansion  and  free 
bubbling  begins  at  or  a  little  above  minimum  fluidization.  Geldart  states 
that  a  particle  would  generally  belong  to  Group  B  if 

(Ps  -  0g)1J7  dp  >  906,000  (2.4 

with  the  densities  expressed  in  kg/m^  and  the  particle  diameter  in  ym. 
Group  C  materials  exhibit  a  strong  tendency  to  agglomerate  and  are  usually 
difficult  to  fluidize.  Particles  in  Group  D  generally  fluidize  unstably. 
In  air  at  atmospheric  pressure  a  particle  would  belong  to  Group  D  if 

(ps  -  Pg)dp2  >  109  (2.5 

2.2  The  Mechanisms  of  Heat  Transfer  in  a  Fluidized  Bed 

In  an  aggregatively  fluidized  bed,  it  is  generally  believed  heat 
transfer  is  caused  by  the  rapid  mixing  of  the  particles  within  the  bed 
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FIGURE  2.1  Geldart's  Materials  Classification  [5] 
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Bubbles  displace  particles  as  they  pass  and  entrain  particles  into  their 
trailing  wake.  As  a  bubble  sweeps  near  or  over  a  surface,  particles  at  the 
surface  are  displaced  by  fresh  particles  at  the  mean  temperature  of  the 
bed.  Heat  transfer  ensues  as  long  as  there  is  a  temperature  difference 
between  the  particles  and  the  surface.  It  has  been  suggested  that 
initially  the  heat  transfer  is  confined  to  the  surface  and  first  layer,  of 
particles  adjacent  to  it  [6].  As  time  progresses,  the  temperature  of  the 
adjacent  particles  approaches  that  of  the  surface  and  heat  transfer  ensues 
between  the  first  and  second  layer  of  particles.  If  the  particles  remain 
at  the  wall  for  a  long  period  of  time  the  heat  transfer  interaction 
includes  particles  far  from  the  wall. 

Mickley  and  Fairbanks  [7]  were  the  first  to  physically  model  a 
transient  renewal  model  of  heat  transfer  between  an  immersed  surface  and  a 
fluidized  bed.  In  this  model  the  dense  phase  (the  mixture  of  particles  and 
gas)  adjacent  to  the  surface  is  considered  a  homogeneous  material  with  an 
effective  thermal  conductivity  and  heat  capacity.  Mickley  and  Fairbanks 
postulated  that  heat  transfer  ensues  as  a  packet  of  this  material  comes  in 
contact  with  the  immersed  surface.  A  one-dimensional  transient  heat 
conduction  takes  place  until  the  packet  is  replaced  by  another  one  due  to 
bubble  motion  sweeping  the  surface  or  otherwise  disrupting  the  material. 

The  resulting  heat  transfer  coefficient,  hg  ,  predicted  by  this  model 
averaged  over  the  packet  residence  time  at  the  surface  is 
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where 

kg  =  effective  packet  thermal  condutivity 

Ps(l-em^)Cs  =  effective  packet  heat  capacity 

tr  =  packet  residence  time  at  the  surface. 

This  model  of  heat  transfer  predicts  that  as  the  residence  time 
continually  decreases  heat  transfer  will  increase.  However,  several 
investigators  have  experimentally  found  that  the  heat  transfer  coefficient 
approaches  an  upper  limit  as  residence  times  are  decreased  [8,9].  Baskakov 
[10]  accounted  for  this  result  by  introducing  the  concept  of  a  thermal 
constant  resistance  in  series  between  the  surface  and  particles. 

Additional  work  has  further  indicated  that,  especially  for  large 
particles  (~lmm),  heat  transfer  between  particles  and  a  surface  is  not  only 
dependent  upon  the  heat  capacity  of  the  particle  and  the  direct  thermal 
contact  between  particle  and  surface,  but  also  on  convective  and  radiative 
heat  transfer  through  the  gas  layer  surrounding  the  contact  points  [6]. 
Further  it  has  been  demonstrated  that  the  average  heat  transfer  coefficient 
for  large  particle  beds  shows  a  modest  rise  with  increasing  U/Um^ 

approaching  a  constant  value.  Small  particle  beds,  on  the  other  hand,  have 
exhibited  dramatic  changes  in  heat  transfer  coefficient  as  U/Umf 

increases  [6], 

In  light  of  these  results,  Gelperin  and  Einstein  [11]  have  postulated 
that  the  heat  transfer  averaged  over  an  immersed  surface,  h  ,  is 
determined  by  the  percent  of  the  surface  covered  by  particles,  (1-6),  the 
average  heat  transfer  coefficient  at  the  surface  covered  by  the  emulsion, 
^surface’  an<*  avera9e  lleat  transfer  coefficient  during  bubble  or  void 
contact  at  the  surface,  h  , 
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In  this  model  the  average  heat  transfer  between  the  particles  and  the 
surface,  hsurface  »  1S  due  to  the  previously  discussed  combined  series 
effect  of  the  surface  resistance  to  heat  transfer  and  the  particle  renewal 
model  of  Mickley 
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where 


=  the  combined  effect  of  conduction,  convection 
and  radiation  heat  transfer  at  the  surface 
=  the  heat  transfer  coefficient  defined  by  Eq. 
(2.6). 


In  non-dimensional  terms  h  can  be  expressed  as 
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and  using  values  determined  by  Glicksman,  et  al.  [6] 
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The  value  of  6  in  Eq.  (2.10a)  is  based  on  Decker's  [1]  results  in  a  well 
fluidized  bed.  The  value  of  12  in  Eq.  (2.10b)  is  based  upon  Gloski's  [12] 
results  for  near  a  bed  wall  at  close  to  minimum  fluidization  conditions. 

Although  the  bubbles  in  a  fluidized  bed  help  promote  heat  transfer  by 
providing  the  mechanism  for  rapid  particle  replacement  against  an  immersed 
surface,  they  also  tend  to  reduce  the  overall  heat  transfer  by  covering  a 
portion  of  the  surface  and  thus  reducing  the  surface  area  in  contact  with 
particles.  This  effect  is  represented  by  the  term  h  ..  in  Eq.  (2.7). 

The  non-dimensional  expression  for  may  be  estimated  using  the 

Polhausen  solution  for  heat  transfer  in  a  laminar  boundary  layer  on  a  flate 
plate 


Nu  .  ,  =  0.664  Re0,5  Pr0,33  (  i  )°’5  (2.11) 
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where  Re  =  the  Reynolds  Number  expressed  as  [6] 
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Pr  =  the  gas  Prandtl  Number 

d^  =  the  average  bubble  diameter  passing 


over  the  surface. 


The  model  represented  by  Eq.  (2.7)  has  been  used  by  several 
investigators  [1,6,12]  to  predict  heat  transfer  coefficients  with  good 
results  and  will  be  used  in  this  investigation. 

Another  area  of  potentially  affecting  fluidized  bed  heat  transfer  is 
through  the  use  of  finned  surfaces.  The  presence  of  fins  increases  the 
available  surface  area  for  heat  transfer.  This  may  be  especially  useful 
for  water  cooled  immersed  tubes  when  the  major  resistance  to  heat  transfer 
occurs  on  the  outside  of  the  tube. 

However,  there  presently  appears  to  be  a  general  lack  of  extensive 
data  available  for  finned  surface  use  in  a  fluidized  bed.  This  lack  of 
data  precludes  the  use  of  any  available  model  or  correlation  to  predict 
finned  surface  heat  transfer  coefficients  with  confidence.  Glicksman  has 
suggested  [6]  using  data  for  bare  tubes  to  predict  h  of  Eq.  (2.7)  and 
then  reducing  this  value  by  30-50%  to  estimate  the  finned  tube  heat 
transfer  coefficient.  Further,  it  has  been  shown  by  various  investigators 
[13,14,15,16]  that  in  spite  of  this  reduction  in  heat  transfer  coeficient 
and  the  inclusion  of  fin  efficiency,  finned  surfaces  still  can  provide 
increased  overall  heat  transfer.  It  has  also  been  found,  though,  that  as 
the  fin  spacing  is  reduced  to  less  than  ten  particle  diameters  the  heat 
transfer  coefficient  decreases  due  to  the  probable  decreased  particle 
mobility  between  the  fins  [6]. 


CHAPTER  3 


CONCEPT  OF  THE  INVESTIGATION 


Chapter  2  outlined  the  general  theory  and  relationships  that  have  been 
developed  to  study  particle  fluidization  and  fluidized  bed  heat  transfer  in 
air.  This  chapter  will,  initially,  define  the  specffic  areas  of  interest 
of  this  investigation.  A  further  development  of  the  concepts  presented  in 
Chapter  2  will  then  follow  in  order  that  they  may  be  directly  applicable  to 
the  requirements  of  this  work. 

3.1  Areas  of  Investigation 

Chapter  1  discussed  the  need  for  analyzing  the  fluidization  and  heat 
transfer  characteristics  of  low  density  particles  in  a  fluidized  bed.  In 
this  investigation  several  particles  with  densities  of  approximately  2-30 
pcf  and  diameters  of  0.25  -  2  mm  will  be  analyzed  by  attempting  to  apply 
the  theory  discussed  in  Chapter  2.  Experiments  will  be  conducted  in  which 
fluidization  and  heat  transfer  data  will  be  collected  and  analyzed. 

Fluidization  data  of  interest  in  this  investigation  is  the  minimum 
fluidization  velocity  of  the  test  particles,  the  pressure  drops  associated 
with  the  low  density  aggregatively  fluidized  beds,  and  the  Geldart 
classification  of  the  test  particles  based  upon  their  ability  to  be 
fluidized.  Heat  transfer  data,  in  the  form  of  average  heat  transfer 
coefficients,  will  be  collected  for  two  geometries:  wall-to-bed  and 
immersed  tube-to-bed.  Relationships  between  heat  transfer  coefficients  and 
air  velocity,  particle  density,  and  particle  diameter  will  be  investigated. 
Comparisons  between  unobstructed  and  obstructed  air  flow  on  wall  heat 
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transfer  and  between  finned  and  unfinned  immersed  tube-to-bed  heat  transfer 
are  also  of  interest. 

Although  the  present  theory  available  for  this  investigation  has  been 
developed  using  significantly  denser  material  it  will,  nevertheless,  be 
utilized  as  a  basis  for  comparing  results  of  low  density  material  testing. 
As  such,  the  fluidization  theory  outlined  in  Chapter  2  can  be  compared 
directly.  The  heat  transfer  theory  will  require  further  detailing  to  be 
satisfactorily  used  for  this  analysis.  Section  3.2  presents  this  detail. 
3.2  Model  for  Heat  Transfer  Analysis 

The  model  for  heat  transfer  analysis  that  is  used  in  this 
investigation  is  the  one  expressed  by  Eq.  (2.7).  Inherent  in  this  model 
are  several  factors  that  must  be  evaluated  or  reasonably  assumed  based  upon 
the  model's  particular  application. 

The  first  factor  entails  the  determination  of  the  fraction  of  a 
surface  covered  by  voids,  6  .  There  is  no  specifically  accurate  way  of 
determining  6  .  However,  several  methods  and  relationships  exist  for 
estimating  it  closely.  This  analysis  will  be  based  on  Geperin's  [11] 
relationship  of 


6  =  1  -  ( 
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with  e  and  e  ^  having  assumed  values  of  0.6  and  0.5,  respectively,  for 
all  test  particles.  These  values  are  based  upon  representative  data  from 
other  investigations  of  comparable  particle  diameter  but  greater  density 
fluidization  [17].  Substituting  these  values  for  e  and  into  Eq. 

(3.1)  yields 
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The  term  representing  the  average  heat  transfer  between  the  particles 
and  the  surface,  hsurface  »  is  defined  by  Eq.  (2.8).  Contained  within 
this  expression  are  several  parameters  represented  by  Eqs.  (2.6)  and 
(2.10).  The  heat  transfer  coefficient  of  the  emulsion,  hfi  ,  theorized  by 
Mickley  is  based  upon  effective  thermal  properties  and  parti  cl e-to-surf ace 
residence  times.  These  parameters  are  inherently  difficult  to  quantify  and 
are  themselves  subject  of  extensive  research  [17].  In  this  investigation 
Gelperin's  relation  for  effective  conductivity,  kg  ,  will  be  used  [11]: 
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Parti  cl e-to-surface  residence  times  will  be  estimated  based  upon  visual 
observation  and  previous  experimental  representative  values  [17].  This 
inherent  uncertainty  of  residence  time  is  a  possible  critical  area  of 
analysis  and  will  be  further  discussed  in  Chapter  6. 

The  heat  transfer  coefficient  represented  by  h^  can  be  derived  from 
Eq.  (2.10) 
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where  Re  =  Reynolds  Number  defined  as  [6] 
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The  radiation  component  Nusselt  Number,  Nupa(j  ,  may  be  expressed  by  [6] 
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=  emissivity  of  the  fluidized  bed 
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Combining  Eqs.  (3.4a)  and  (3.4b)  with  Eq.  (3.5)  and  substituting  the 
appropriate  air  property  values  the  following  relations  for  h,.  result 
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Both  expressions  for  h^  will  be  used  in  this  investigation  for  analysis 
of  the  low  density  particle  heat  transfer  due  to  the  unknown  nature  of 
their  results. 

The  term  representing  the  heat  transfer  due  to  the  passage  of  bubbles 
or  voids  over  a  surface,  h  »  1S  defined  in  nondimsional  terms  by  Eq. 
(2.11).  Further  quantification  of  this  expression  can  be  accomplished 
through  substitution  of  known  and  approximated  parameters.  The  average 
bubble  rise  velocity,  ,  is  closely  approximated  by  [6] 


Ub  =  0.71^  ♦  (U0  -  Umf) 


(3.7) 


For  the  purpose  of  this  analysis  and  to  be  conservative  the  average  bubble 
length,  db  ,  will  be  taken  as:  one-half  the  plate  width  for  the  wall 
heater,  one-quarter  of  the  tube  circumference  for  the  unfinned  tube  heater, 
and  equal  to  the  average  fin  spacing  for  the  finned  tube  heater.  Using 
these  approximations  and  the  appropriate  air  properties  h  ^  may  be 
expressed  as 


Un  +  2U  ,  +  0.71/g~3T 

u  _  „  ,n,r  0  mf  3  b  n 
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Finally,  for  finned  surface  analysis  two  approaches  will  be 
investigated.  The  first  approach  will  be  to  simply  apply  the  previously 
described  heat  transfer  model  modified  by  an  overall  corresponding  surface 
efficiency.  For  a  finned  surface  the  total  rate  of  heat  transfer,  q  ,  is 


Taking  hu  =  hf  =  ^  *  an  overan  surface  efficiency,  e  ,  may  be  defined 

as 
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where  A  =  Au  +  Af  [18]. 

The  second  approach  will  be  to  numerically  compare  the  finned  tube 


heat  transfer  coefficients  to  the  unfinned  values  to  determine  the 
magnitude  of  change,  if  any,  between  the  two  geometries. 
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CHAPTER  4 

EXPERIMENTAL  APPARATUS  DESIGN 

Although  the  proposed  fluidization  and  heat  transfer  experiments  are 
straightforward  in  theory,  the  design  and  construction  of  the  various 
components  of  the  experimental  apparatus  requires  care.  This  chapter  will 
discuss  these  aspects  of  the  experimental  apparatus  employed  to  carry  out 
the  necessary  data  collection.  This  discussion  will  include  the  fluidized 

beds,  heaters,  and  particles  used  in  these  experiments.  The  details  of  the 

data  acquisition,  however,  will  be  discussed  in  Chapter  5. 

4.1  Fluidized  Beds 

Two  different  fluidized  beds  have  been  used  for  data  collection.  A 
two-dimensional  bed  is  used  primarily  to  study  the  fluidization  behavior  of 
the  test  particles  and  the  wall  heat  transfer  characteristics.  This  type 
of  bed  allows  experiments  to  be  carried  out  in  a  fluidized  bed  of  similar 
geometry  to  that  of  a  typical  building  wall  cavity.  It  consits  of  a  1  in  x 
16  in  rectangular  cross-section  and  an  overall  height  of  24  inches 
available  for  particle  fluidization.  The  air  distributor  is  a  plastic 
perforated  sheet  covered  underneath  by  three  inches  of  foam  rubber.  The 
foam  rubber  serves  to  increase  the  pressure  drop  across  the  distributor  and 

thus  insure  a  more  uniform  air  flow.  Figure  4.1  details  the  overall 

dimensions  of  the  two-dimensional  bed. 

A  three-dimensional  bed  is  used  to  study  the  immersed  tube  heat 
transfer  characteristics.  Although  this  type  of  bed  does  not  realistically 
model  the  rectangular  cross  section  of  a  typical  wall  cavity,  this 
particular  bed  allows  for  good  particle  to  tube  interaction  that  is  not 


possible  in  the  available  two-dimensional  bed.  The  depth  of  the 
three-dimensional  bed  and  the  radial  dimension  of  the  immersed  tube  heaters 
approximates  that  of  the  depth  of  a  wall  cavity,  3-4  inches,  and  a  standard 
immersed  heat  transfer  tube,  1/2  to  3/4  inch  in  diameter,  within  it.  In 
this  arrangement  there  is  approximately  a  2-1/2  inch  clearance  between  the 
unfinned  tube  heater  and  the  bed  walls  (front  and  back)  and  a  1-7/8  inch 
clearance  between  the  finned  tube  heater  and  the  walls.  This  bed  consists 
of  a  6  in.  x  6  in.  square  cross-section  and  an  overall  height  of  36  inches 
available  for  particle  fluidization.  The  air  distributor  consists  of  a 
porous  plate  covered  by  approximately  16  inches  of  sand  to  increase  the 
pressure  drop  across  the  distributor  plate  and  insure  uniform  air  flow 
through  the  bed.  Figure  4.2  details  the  overall  dimensions  of  the 
three-dimensional  bed. 

The  flow  of  air  through  both  beds  is  metered  by  a  standard  ASME 
square-edged  orifice  with  flange  taps  [19].  Each  bed  has  a  tap  near  the 
bottom  for  particle  removal  and  a  plexiglas  front  wall  for  ease  of  particle 
observation.  Additionally,  both  beds  have  pressure  taps  located  on  the 
walls  available  for  measuring  pressure  drops  within  the  bed. 

4.2  Flat  Plate  (Wall)  Heater  Design 

The  measurement  of  average  heat  transfer  coefficients  from  a  wall 
surface  to  the  fluidized  ma trial  is  accomplished  by  the  use  of  a  flat  plate 
wall  heater.  To  insure  proper  bed  material  fluidization  and  heat  transfer 
at  the  bed  walls  it  is  necessary  for  the  wall  heater  to  be  designed  such 
that  the  heat  flux  generated  by  the  heater  is  transferred  directly  into  the 
two-dimensional  bed  and  not  lost  to  other  areas,  the  heater  surface  remains 
at  a  uniform  temperature,  and  the  heater  is  mounted  flush  with  the  existing 
bed  wall  thus  maintaining  a  smooth  profile. 


Since  the  objective  of  the  experiments  using  the  wall  heater  is  to 
measure  the  heat  transfer  coefficient  from  a  heated  wall  surface  to  a 
fluidized  bed  of  low  density  particles,  it  is  vital  that  the  experimenter 
know  what  heat  flux  is  going  into  the  bed.  Thus,  heat  losses  from  other 
areas  of  the  heater  not  in  contact  with  the  fluidized  particles  are  to  be 
kept  to  a  minimum.  Further,  a  uniform  wall  temperature  must  also  be 
maintained  on  the  heater  surface  in  order  to  accurately  calculate  the  heat 
transfer  coefficient,  and  the  heater  element  must  be  positioned  in  the 
bed's  existing  wall  in  order  to  both  properly  transfer  the  desired  heat 
flux  and  not  affect  the  normal  fluidizing  characteristics  of  the  test 
material.  Unwanted  protrusions  from  the  wall  into  the  bed  caused  by  the 
heater  may  decrease  the  bed  materials  mixing  and  consequently  lower  the 
heat  transfer  coefficient. 

A  flat,  well  insulated  heater  made  of  high  thermal  conductivity 

material  is  used  in  this  experiment.  As  seen  in  Figure  4.1,  the  front 

plexiglas  wall  of  the  two-dimensional  bed  measures  16  in  x  24  in.  Due  to 

the  size  of  this  frontal  area  a  design  choice  has  been  made  to  confine  the 

heater  size  to  1.5  in  x  6  in.  This  size  provides  an  adequate  heat  transfer 
2 

surface  (0.063  ft  )  while  still  allowing  visual  observation  of  the  material 
being  fluidized  within  the  bed.  The  heater  is  oriented  in  the  plexiglas 
wall  with  the  long  dimension  being  vertical.  This  allows  for  the  longest 
heating  length  to  be  in  the  same  direction  as  the  bed's  air  flow  and  bubble 
direction  of  travel. 

Secondly,  a  0.25  in.  thick  copper  plate  was  chosen  as  the  heated 
surface  due  to  copper's  high  thermal  conductivity  and  to  insure  good 
temperature  uniformity  on  the  heater  surface  in  contact  with  the  bed 
material.  Styrofoam  is  used  as  the  primary  heater  insulation  due  to  its 


good  insulating  properties  and  availability.  A  standard  150  watt,  120  volt 
strip  heater  is  used  as  the  heat  source  to  be  attached  to  the  copper  plates 
by  two  0.25  in  x  1  in  threaded  brass  rods.  Holes  are  drilled  in  the  copper 
plate  approximately  0.0313  in.  in  from  the  particle  contacting  surface  for 
copper/constantan  thermocouple  placement. 

Dimensions  for  the  insulation  required  to  insure  small  heat  losses, 

on  the  order  of  1%  of  the  total  heat  flux,  were  determined  through  a  heat 

balance  around  the  heater.  For  these  calibrations  an  assumed  value  of  the 

heat  transfer  coefficient  from  the  heated  surface  to  the  bed  of  20 
2 

Btu/hr-ft  -°F  and  the  ratio  of  the  difference  between  the  wall  and 
fluidized  bed  temperatures  to  the  difference  between  the  wall  and  ambient 
air  temperatures  of  0.20  is  used.  As  will  be  seen  in  Chapter  5  these 
assumptions  are  slightly  over  optimistic  but  do  not  compromise  the  results 
of  the  design.  Appendix  A  outlines  the  procedure  used  in  determining  the 
insulation  dimensions  via  the  heat  balance  analysis  and  Figure  4.3  depicts 
the  final  design  of  the  wall  heater. 

4.3  Unfinned  Immersed  Tube  Heater  Design 

The  experimental  requirements  for  this  heater  are  similar  to  those  of 
the  wall  heater.  However  its  design  and  construction  are  somewhat  easier 
due  to  the  drastic  reduction  of  insulation  required.  The  purpose  of  this 
heater  is  to  provide  a  known  heat  flux  into  the  three-dimensional  fluidized 
bed  and  maintain  a  uniform  tube  wall  temperature  while  being  immersed 
within  the  bed. 

For  this  design  caution  should  be  exercised  in  order  to  insure  good 
thermal  contact  between  the  heating  element  and  the  tube  material. 

Although,  with  proper  insulation  heat  losses  will  be  minimal,  poor  thermal 
contact  between  the  heating  element  and  tube  and/or  thin  tube  wall 
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FIGURE  4.3  Wall  Heater  Design 


thickness  could  lead  to  nonuniform  tube  surface  temperatures.  This,  in 
turn,  will  make  accurate  calculation  of  heat  transfer  coefficients 
difficult.  For  this  reason  standard  wall  tubing  is  not  used  in  this 
experiment.  Rather,  a  1  in  (diameter)  x  5  in  aluminum  bar  has  been 
fabricated  into  the  required  tube  heater. 

A  hole  of  0.375  in  is  drilled  and  reamed  through  the  center  of  the 
aluminum  bar.  Two  additional  holes  on  opposite  ends  of  the  bar  are  drilled 
approximately  0.0313  in  in  from  the  particle  contacting  surface  for 
copper/constantan  thermocouple  placement.  A  0.375  in  x  5  in,  280  watt,  120 
volt  standard  cartridge  heater  is  inserted  into  the  center  hole  of  the 
aluminum  bar.  Styrofoam  insulation  is  placed  on  the  ends  of  the  bar  to 
prevent  heat  loss  and  provide  structural  stability  for  placement  in  the 
bed. 

This  design  utilizes  a  high  thermally  conductive  material  and  provides 
for  over  0.25  inches  of  tube  wall  thickness  to  insure  uniform  tube  surface 
temperature.  Its  overall  length  allows  it  to  be  centered  and  rigidly 
placed  in  the  three-dimensional  bed  thus  preventing  any  possible  adverse 
effects  to  the  bed  material  fluidization  and  heat  transfer  due  to  heater 
motion  during  fluidization.  Figure  4.4  depicts  the  design  of  the  unfinned 
immersed  tube  heater. 

4,4  Finned  Immersed  Tube  Heater  Design 

This  heater  only  differs  from  the  unfinned  tube  heater  in  concept  by 
the  presence  of  fins.  Once  again,  however,  standard  finned  tubing  presents 
the  same  thin  tube  wall  and  heating-element- to- tube- thermal  contact 
problems  that  are  discussed  with  the  unfinned  tube  design.  Additionally, 
in  many  cases  standard  finned  tubing  has  questionable  to  poor  thermal 
contact  between  the  fin  base  and  the  tube.  It  may  be  that  in  actual 


application  standard  finned  tubing  will  be  used  (see  Chapter  7),  but  the 
purpose  of  this  heater  is  to  accurately  assist  in  the  determination  of  a 
heat  transfer  coefficient.  It  is  therefore  necessary  to  reduce  all 
possible  sources  of  error  and  misrepresentation  by  insuring  good  thermal 
contacts. 

For  these  reasons,  the  choice  has  been  made  to  fabricate  a  finned 
tube  heater.  The  design  of  this  unit  consists  of  a  0.75  in  x  4.5  in  copper 
bar  with  13  copper  circumferential  fins  of  rectangular  cross-section 
attached.  A  0.25  in  hole  is  drilled  and  reamed  through  the  center  of  the 
bar.  Thermocouple  holes  are  placed  similarly  to  those  on  the  unfinned 
tube.  A  0.25  in  x  4.5  in,  110  watt,  120  volt  standard  cartridge  heater  is 
inserted  into  the  center  hole  of  the  copper  bar.  Styrofoam  insulation  is 
placed  on  the  ends  of  the  bar  to  prevent  heat  loss  and  provide  stability  in 
the  bed. 

As  with  the  unfinned  tube,  a  high  thermally  conductive  material  with 
approximately  0.25  in  tube  wall  thickness  is  used.  A  copper  tube  with 
copper  fins  is  used  rather  than  aluminum,  as  with  the  unfinned  tube, 
because  of  the  need  to  provide  for  good  fin  to  tube  thermal  contact. 

Copper  on  copper  can  be  easily  and  thermally  effectively  soldered  together, 
whereas  aluminum  cannot. 

The  copper  fins  are  2.25  inches  in  diameter,  square  edged,  and  0.02 
inches  thick.  They  are  spaced  approximately  3  fins  per  inch  on  the  tube 
and  have  an  efficiency  of  90%.  The  fins  are  soldered  to  the  copper  tube 
with  standrd  50/50  (Pb-Sn)  solder.  Care  must  be  taken  to  achieve  good 
thermal  contact  around  the  fin  base  with  the  tube.  Appendix  B  details  the 
specific  analysis  of  these  fins.  Figure  4.5  illustrates  the  finned 
immersed  tube  design. 


4.5  Particle  Description 

The  primary  distinction  and  sole  purpose  of  this  investigation  is  to 
analyze  the  fluidization  and  heat  transfer  of  low  density  particles  in  a 
fluidized  bed.  As  was  discussed  in  Chapter  2,  previous  work  has  been 
extensively  done  on  high  density  material  in  a  fluidized  bed  and  the 
results  of  their  fluidization  and  heat  transfer  have  been  well  documented 
[12,20,21]. 

Several  different  particles  have  been  tested  in  this  experiment  with 
various  results.  These  results  are  discussed  in  Chapter  5  and  6.  Seven 
specific  particles  are  thoroughly  investigated.  They  range  in  average 
particle  diameter  of  250  ym  to  2  mm  and  density  from  approximately  2.5  pcf 
to  32  pcf. 

Test  particles  used  have  come  from  three  primary  sources:  the  3M 
Company  [22],  the  Norton  Company  [23],  and  the  Whittermore  Perlite  Company 
[24].  The  3M  "Macrospheres"  are  hollow  cells  with  a  shell  made  of  a 
hollow  glass  bubble/resin  combination.  Their  present  applications  include 
fiberglas  reinforced  plastics  construction,  filters,  flotation  devices, 
marine  fabrication,  vapor  suppression,  and  electrolytic  spacers.  In  this 
application  the  3M  "Macrospheres"  M40X,  particle  diameter  of  1  mm,  density 
of  27.2  pcf  and  the  M27X,  particle  diameter  of  2  mm,  density  of  15.31  pcf 
are  used. 

Three  particles  are  used  from  the  Norton  Company.  Thv;se  particles  are 
Norton  650,  diameter  of  650  ym,  density  of  31.78  pcf,  Norton  1,  diameter  of 
Irnrn,  density  of  23.2  pcf,  and  Norton  1.5,  diameter  of  1.5mm,  density  of 
20.85  pcf.  These  cellular  glass  nodules  are  composed  of  hundreds  of 
individual  hermetically  sealed  closed  cells.  Their  present  applications 
include  use  as  a  lightweight  insulating  filler  in  composites  and  for 
inducing  controlled  pores  in  ceramic  bodies. 
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The  third  group  of  test  particles  are  two  sizes  of  Perlite.  Perlite 
is  a  generic  term  for  a  naturally  occurring  siliceous  volcanic  rock.  When 
heated  to  a  suitable  point  in  its  softening  range,  it  expands  four  to 
twenty  times  its  original  volume.  This  expansion  is  due  to  the  presence  of 
two  to  six  percent  combined  water  in  the  crude  perlite  rock.  Upon  heating, 
the  water  vaporizes  creating  microscopic  hollow  glassy  bubbles.  The 
formation  of  these  bubbles  in  the  perlite  causes  the  expansion  and  accounts 
for  perlite's  characteristic  low  density.  The  expansion  process  can  be 
controlled  allowing  for  the  manufacturing  of  perlite  ranging  in  bulk 
density  from  2-25  pcf  and  particle  diameter  of  250  pm  to  4mm.  Its  present 
applications  are  strongly  in  use  in  the  building  industry  as  a  thermal 
insulator,  acoustical  plastic  ingredient,  and  a  concrete  aggregate  and  in 
the  agricultural  industry  as  a  soil  conditioner  and  fertilizer  extender. 

In  the  present  application  of  perlite  in  this  investigation  the  average 
particle  diameter  is  250  pm,  with  a  density  of  4.5  pcf,  and  a  diameter  of 
500  m,  and  a  density  of  2.7  pcf. 

For  all  particles  in  this  investigation,  the  stated  particle  diameter 
represents  the  average  diameter  determined  from  a  standard  seive  analysis 
and  may  vary  by  +  20%.  Densities  are  expressed  as  bulk  particle  density. 
They  have  been  determined  by  first  measuring  the  loose-fill  density  of  the 
particles  and  then  dividing  this  value  by  the  percent  volume  occupied  by 
the  solid  particles.  For  this  calculation  all  the  particles  were  assumed 
to  be  spherical  in  shape  and  a  loose  fill  voidage  of  0.47  was  used.  Thus 
bulk  density  was  thus  determined  by 

Loose-fill  Density 


Bulk  Density 


(4.1) 


CHAPTER  5 


DATA  ACQUISITION  AND  REDUCTION 

In  the  last  chapter  the  design  and  construction  of  the  required 
experimental  aparatus  and  a  description  of  the  low  density  particles  used 
was  conducted.  This  chapter  will  detail  the  procedures  used  in  acquiring 
the  fluidization  and  heat  transfer  data  of  the  seven  test  materials,  the 
reduction  of  this  data  to  usable  form  for  later  analysis,  the  checks  made 
to  determine  heater  performance,  and  an  estimate  of  the  error  associated 
with  the  acquired  data. 

5.1  Particle  Fluidization  Determination 

Since  no  previous  work  has  been  conducted  on  the  actual  fluidization 
of  low  density  material,  this  determination  has  to  be  made  before  heat 
transfer  experiments  can  be  performed.  As  described  in  Chapter  4,  seven 
particles  have  been  completely  analyzed  for  this  purpose.  However,  there 
have  been,  at  least  initially,  more  than  these  seven  particles 
investigated.  In  this  section  the  inappropriateness  of  these  other 
particles  and  the  determination  of  the  fluidization  characteri sties  of  the 
seven  test  particles  is  discussed. 

To  qualify  as  a  test  particle,  the  particle  must  exhibit  good 
fluidization  characteristics.  It  must  be  capable  of  being  fluidized  as 
evidenced  by  the  presence  of  good  bubble  motion,  good  particle  mixing,  and 
modest  air  pumping  power  requirements  for  fluidization.  The  seven 
previously  described  materials  satisfy  these  requirements,  but  not  all 
particles  investigated  did. 

Three  such  additional  categories  of  particles  have  been  identified. 

The  first  category  is  of  the  group  that  is  characterized  by  a  large  average 


particle  diameter  and  small  density.  A  material  illustrating  this  is  0.125 
in  (3  mm)  styrofoam  beads  having  a  density  of  approximately  2.35  pcf.  When 
placed  in  the  two-dimensional  bed  it  is  visually  observed  that  fluidization 
is  very  poor.  The  beads  appear  to  exhibit  an  electrostatic  attraction  to 
the  walls  of  the  bed.  There  appears  to  be  no  noticeable  attraction  between 
the  beads  themselves  and  some  fluidization  takes  place  in  the  interior  of 
the  bed  away  from  the  walls.  '  Due  to  the  insulating  nature  of  styrofoam  and 
the  static  electric  attraction  to  walls,  this  material  is  not  appropriate 
for  the  intended  experiment.  This  is  not  to  say,  however,  that  these 
particles,  nor  styrofoam  in  general,  should  be  totally  ignored.  As 
discussed  in  Appendix  C  static  charging  is  a  function  of  certain  material 
properties  and  both  styrofoam  and  plexiglas  exhibit  a  strong  electrostatic 
charging  potential.  Procedures  can  be  taken  to  help  reduce  the  effect  of 
this  phenomenon,  a  few  of  which  are  discussed  in  Chapter  8.  These 
procedures  are  not  employed  in  this  experiment. 

Another  category  investigated  is  that  of  large  diameter  and  moderate 
density  particles.  As  described  in  Chapter  2,  a  material's  minimum 
fluidization  velocity  is,  among  other  things,  a  function  of  particle 
diameter  and  density.  As  illustrated  by  the  results,  shown  in  Chapter  6, 
for  the  moderate  to  low  density  material  the  minimum  fluidization  velocity 
tends  to  be  a  much  stronger  function  of  particle  diameter  than  density  at 
moderate  densities.  Consequently,  larger  particles  will  generally  require 
more  air  velocity  to  reach  fluidization.  This  greater  amount  of  air 
velocity  directly  relates  to  an  increased  external  air  pumping  power 
requirement  based  on  particle  diameter  alone.  An  attempt  was  made  to 
fluidize  Macrosphere  and  Norton  particles  of  diameter  2.8  and  3  mm  and 
density  from  14-18  pcf  in  the  available  beds.  The  air  velocity  and 


consequently  the  air  pumping  power  required  for  fluidization  of  this 
material  is  far  beyond  what  is  desireable  for  the  intended  applications  of 
this  work.  On  the  barely  acceptable  end  is  the  2  mm  macrosphere  M27X.  Its 
data  will  help  illustrate  the  problem  with  increasing  particle  diameter. 

The  third  category  of  particles  initially  considered  were  particles 
with  very  small  diameters  and  low  density.  Particles  ranging  in  size  from 
20  -  200  urn  (mean  size  of  70  urn)  and  density  of  11  pcf  have  been  tested. 
This  silicate-based  hollow  microsphere  material  exhibits  very  poor 
fluidization  and  generally  resembles  the  description  of  Geldart's  Group  C 
materials,  but  are  best  described  in  his  Group  A  [5].  This  result  is 
further  discussed  in  Chapter  6. 

Thus,  for  the  reasons  discussed,  these  three  materials  are  not 
completely  investigated  in  this  work.  Their  value,  however,  is  in 
illustrating  classes  of  low  density  materials  that,  for  one  reason  or 
another,  are  not  suitable  for  the  intended  applications  of  this  work.  That 
is  to  say,  they  fail  to  adequately  demonstrate  the  criteria  established  in 
Chapter  1  of  exhibiting  good  fluidization  characteristics,  having  a  low  bed 
pressure  drop,  and  requiring  modest  external  air  pumping  power.  On  the 
other  hand,  the  seven  previously  described  test  particles,  based  on  initial 
screening,  appeared  to  meet  the  established  criteria  and  thus  qualified  for 
further  investigation.  A  description  of  their  investigation  is  now 
presented. 

To  determine  the  fluidization  characteristics  of  the  test  particles 
the  following  procedures  have  been  used.  All  air  velocities  are  determined 
from  the  use  of  an  ASME  square  edged  orifice  with  flange  taps.  Two  orifice 
plates  have  been  used  throughout  the  experiment  to  allow  for  ease  of 
measurement  and  are  described  in  Appendix  D.  In  the  two-dimensional  bed 


the  particles  under  1  mm  in  diameter  have  used  orifice  plate  A  and  the 
larger  particles  have  used  plate  8.  In  the  three-dimensional  bed  plate  B  is 
used  for  all  particles.  Due  to  the  particular  geometry  of  the  fluidized 
beds  and  the  varying  particle  sizes  used  in  this  investigation  this  use  of 
the  two  different  orifice  plates  resulted  in  more  accurate  measurements  of 
the  orifice  plate  pressure  drops. 

All  pressure  drops,  furthermore,  are  measured  in  one  of  two  ways. 

The  pressure  drop  across  an  orifice  plate  is  measured  by  a  water  filled 
U-tube  manometer  graduated  in  tenths  of  an  inch.  Pressure  drops  in  the 
fluidized  beds  are  measured  by  an  inclined  manometer  calibrated  to  read 
hundredths  of  an  inch  of  water.  Static  air  pressure  upstream  of  the 
orifice  is  read  by  a  standard  pressure  gage  in  psi. 

Particle  minimum  fluidization  was  determined  experimentally  using  the 
procedure  outlined  in  Chapter  2.  Well  beyond  minimum  fluidization  the 
pressure  drop  associated  with  the  fluidized  particles  in  the  bed  remains 
relatively  constant  and  independent  of  the  air  velocity  passing  through  the 
bed.  As  the  air  velocity  decreases,  approaching  minimum  fluidization 
velocity  the  bed  pressure  drop  remains  relatively  constant.  Below  minimum 
fluidization,  however,  the  bed  pressure  drop  decreases  linearly  with  air 
velocity.  The  corresponding  air  velocity  associated  with  the  intersection 
of  the  constant  pressure  drop  line  and  the  linearly  decreasing  line  is  the 
minimum  fluidization  velocity.  This  procedure  is  used  to  determine  the 
particle  minimum  fluidization  velocities  and  bed  pressure  drops  of  the  test 
particles.  Figure  5.1  illustrates  this  procedure.  Appendix  E  displays  the 
results  graphically  of  these  determinations  for  the  test  particles. 
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5.2  Determination  of  Heat  Transfer  Coefficients 

To  determine  all  heat  transfer  coefficients,  h  ,  in  this  experiment, 
data  has  been  taken  to  evaluate  or  measure  the  appropriate  terms  contained 
in  the  definition  of  the  heat  transfer  coefficient: 

q  =  h  A  AT  (5-1) 

Measurement  of  all  heat  fluxes  into  the  fluidized  beds  is  accomplished 
by  incorporating  the  heaters  in  a  simple  electric  circuit.  The  voltage  and 
current  supplied  by  a  power  source  to  a  heater  is  measured.  In  this 
experiment  a  variable  A-C  (VARIAC)  power  supply  is  utilized  to  power  the 
heater.  Voltage  and  current  is  measured  independently  by  two  multimeters. 
Using  the  relation 

Power  =  Voltage  x  Current 
or 

q  -  El  (5-2) 

the  heat  flux  may  be  accurately  measured  Figure  5.2  details  the  heater 
electric  circuit. 

The  value  of  each  particular  heater  heat  transfer  area  is  simply  the 
total  surface  area  of  the  heater  in  contact  with  the  fluidized  particles. 

2 

As  indicated  in  Figures  4.3  through  4.5  these  constant  values  are  0.063  ft 

2  2 
for  the  wall  heater,  0.109  ft  for  the  unfinned  tube  heater,  and  0.644  ft 

for  the  finned  tube  heater. 

Temperature  measurements  have  been  taken  using  standard 
copper/constantan  thermocouples  and  read  by  a  digital  copper/constantan 
thermocouple  thermometer  to  a  1°F  accuracy.  Thermocouples  have  been  placed 
in  various  locations  on  the  heaters  and  in  the  fluidized  beds  to  insure 
accurate  temperature  measurement  and  uniformity  on  the  heater  surfaces  and 
in  the  beds  (see  Figures  4. 1-4. 5). 


9 


For  all  experimentation,  wall,  unfinned  and  finned  tube  heating,  a 
minimum  of  three  tests  per  particle  batch  has  been  made.  Heat  flux  for 
each  test  is  kept  approximately  constant  and  the  ratio  of  air  velocity  to 
minimum  fluidization  velocity,  U/U^  ,  is  varied  test  to  test.  Each  test 
is  allowed  to  reach  a  steady-state  condition,  that  is,  the  temperature  of 
the  heater  surface  and  bed,  is  uniform  and  constant  with  time  for  a  given 
heat  flux  and  U/Um.p  before  data  is  taken  and  the  appropriate  heat 
transfer  coefficient  is  evaluated.  To  reach  this  steady  condition  takes, 
on  the  average,  two  hours.  Once  a  steady-state  condition  has  been  reached, 
the  appropriate  data  is  taken  and  the  heat  transfer  coefficient  is 
determined  using  Eq.  (5.1). 

5.3  Heater  Performance  Checks 

In  Chapter  4  and  Appendix  A  details  of  the  heaters  used  in  this 
investigation  and  the  analysis  used  to  minimize  heat  losses  have  been 
discussed.  It  is  further  necessary  to  check  the  heater's  actual 
performance  in  order  to  gain  an  estimate  of  any  error  present  in 
determining  heat  flux  or  any  required  calibration  to  be  done.  This  section 
will  discuss  the  test  performed  on  the  wall  heater  and  tube  heater  to  check 
their  operation. 

The  performance  of  the  flat  plate  wall  heater  was  checked  using  two 
methods.  The  first  was  by  conducting  a  natural  convection  test  of  the 
heater.  The  heater  unit  (the  heater  and  associated  insulation)  was 
oriented  as  depicted  in  Figure  4.1.  A  heat  flux,  q  ,  was  generated  using 
the  experimental  equipment  earlier  described  in  this  chapter.  The 
steady-state  heater  wall  temperature  and  the  ambient  air  temperature  were 
measured.  From  this  data  an  experimental  heat  transfer  coefficient,  h 

of  1.20  Btu/hr-ft  -°F  was  calculated  using  Eq.  (5.1).  This  value  was 
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compared  to  a  theoretical  value  based  upon  the  same  test  conditions.  The 
theoretical  value  is  expressed  as 


htheor.  =  hNC  +  hRAD 


where 

hNC 


hRAD 


=  the  theoretical  natural  convection 
heat  transfer  coefficient  for  a 
vertical  plate  in  air 
=  0.29( AT/L)  [18] 


the  theoretical  radiation  heat  transfer 


coefficient 
r3 


4F  a  r 


mean 


[18] 


(5.3) 


Using  this  expression  with  the  same  measured  parameters  of  the  test  and  an 
assumed  emissivity  of  oxidized  copper  of  0.25  [18]  yields  an  htheor  of 

p 

1.21  Btu/hr-ft  -°F.  This  compares  very  well  with  the  experimentally 
determined  value  having  an  error  of  approximately  1%.  However,  caution 
must  be  exercised  in  relying  heavily  on  the  results  of  this  particular 
check  due  to  the  uncertainty  of  the  assumed  emissivity  value  of  copper. 

The  second  method  utilized  to  check  the  performance  of  the  wall  heater 
consisted  of  orienting  the  heater  unit  in  stationary  air  as  before  but 
covering  the  copper  plate  with  a  material  of  known  thermal  conductivity  and 
thickness.  A  heat  flux  was  generated  and  the  temperature  difference  across 
the  covering  material  was  measured.  This  measured  value  of  AT  is  then 
compared  to  a  theoretically  determined  value  from 


Extruded  polystyrene  foam  insulation  was  obtained  having  a  conductivity 
measured  by  the  supplier  using  a  guarded  hot  plate  traceable  to  NBS 
standards  of  0.017  Btu/hr-ft^-°F  at  75°F  and  a  thickness  of  0.0104  ft. 
[25].  The  insulation  was  placed  over  the  heater,  a  heat  flux  through  the 
heater  was  generated,  and  a  temperature  difference  of  40°F  across  the 
insulation  was  measured  using  copper/constantan  thermocouples. 

Substituting  the  appropriate  known  and  measured  values  into  Eq.  (5.4)  a 
ATf^eor  1S  determined.  Again,  good  agreement  exists  between  the 

actual  and  theoretical  values  to  an  error  of  approximately  1%. 

A  second  check  of  this  method  was  conducted  on  air  flowing  at 
approximately  1.5  fps.  This  air  velocity  was  selected  as  an  estimated 
average  test  velocity  in  the  designed  fluidized  bed  experiments.  Under 
these  conditions,  and  using  the  same  heat  flux  as  before,  a  AT  across  the 
insulation  of  39°F  was  measured.  This  results  in  an  error  of  approximately 
2.8%. 

The  tube  heater  performance  was  similarly  checked  in  stationary  and 
flowing  air  using  the  insulation  covering  method  outlined  above.  In  the 
stationary  air  case  a  AT  of  the  insulation  of  83°F  was  measured.  This 
value  compared  to  the  theoretical  AT  of  84.5°F  yields  an  error  of 
approximately  1.8%.  In  the  flowing  air  case  (Uq  =  1.5  fps)  an  error  of  3% 
resul ted. 

As  a  consequence  of  these  tests  it  has  been  concluded  that  the  heaters 
perform  within  an  acceptable  range  of  experimental  accuracy  and  as 
designed.  Furthermore,  it  appeared  that  any  heat  loss  due  to 


two-dimensional  conduction  effects  is  negligible.  It  is  also  important  to 
note  that  the  k/S,  value  of  the  insulation  (1. 5-2.0  Btu/hr-ft  -°F)  used  in 
these  checks  is  approximately  an  order  of  magnitude  less  than  the  heat 
transfer  coefficient  expected  in  the  fluidized  bed  experiments.  It  is  to 
be  expected,  then,  that  for  the  same  heat  flux  generated  as  in  the  checks 
the  heat  transfer  rate  will  be  greater  for  the  bed  and,  thus,  the  resulting 
error  would  decrease.  A  discussion  of  the  range  of  experimental  accuracy 
is  presented  in  the  next  section. 

5.4  Estimation  of  Experimental  Error 

Section  5.3  has  outlined  the  procedures  and  results  associated  with 
checking  the  heating  performance  of  the  wall  and  tube  heaters.  The  results 
from  this  are,  however,  only  part  the  overall  estimation  of  the 
experimental  error.  This  section  will  outline  the  procedure  used  in 
estimating  this  error. 

There  exists  a  certain  amount  of  inherent  error  in  any  measured  data. 
Since  the  heat  transfer  coefficients,  h,  in  this  work  are  determined  by  Eq. 
(5.1),  the  error  in  the  measured  values  of  the  heat  flux  transferred,  q  , 
the  tube  and  flat  plate  area,  A  ,  and  the  wall  and  bed  temperatures,  Tw 

and  Tg  ,  affect  the  error  in  the  calculated  heat  transfer  coefficient. 
Additionally,  since  air  velocities  through  the  fluidized  beds  and  pressure 
readings  are  also  required,  the  error  involved  in  these  measurements  must 
be  accounted  for.  The  method  that  is  used  to  determine  this  experimental 
error  is  outlined  in  reference  [26]. 

The  accuracy  of  the  meters  used  to  measure  the  voltage  and 
amperage  used  in  determining  q  is  1%  and  2 %  respectively  of  full  scale 
deflection.  Since  q  =  voltage  x  amperage  -  leaks  ,  the  error,  e^  ,  in 
the  measured  q  may  be  calculated  by 


For  both  heaters  the  voltage  and  current  error  band  is  eyg^  =  2% 
and  eAMp  =4%.  As  discussed  in  Section  5.3  the  wall  heater  and  tube 

2.8%  and  3%  as  a  worst  case,  respectively. 


heater  have  an  e 


LEAK 


Consequently,  the  error,  eq  ,  in  the  measured  q  is  5.3%  for  the  wall 


heater  and  5.4%  for  the  tube  heater. 

Heater  dimensions  were  measured  by  a  micrometer,  varied  by  0.4%,  and 
thus  have  a  0.4%  error  band  for  the  area. 

Temperatures  were  measured  by  copper/constantan  thermocouples  to  an 
accuracy  of  +  1°F.  This  results  in  an  error  of  1.4°F  for  the  temperature 
difference  between  the  wall  and  the  bed  material.  Assuming  an  average 
value  of  (Tw  -  T0)  of  60°F  for  both  heaters,  this  results  in  a  2.3%  error 
in  the  measurement  of  the  temperature  difference. 

Therefore,  the  net  error  associated  with  the  determination  of  h  in 
this  work  is  the  square  root  of  the  sum  of  the  squares  of  the  errors  of  the 
individual  contributors  to  h  : 


Error  of  h  =  [ (5.3)2  +  (.4)2  +  (2.3)2]1/2  =  5.8% 

(wall) 


Error  of  h  =  [(5.4)2  +  (.4)2  +  (2.3)2]1/2  =  5.9% 


(tube) 


The  error  associated  in  the  measurement  of  air  velocity  can  be 
estimated  using  the  same  procedure  as  that  used  for  the  heat  transfer 


coefficient.  As  shown  in  Appendix  D,  air  velocity,  Uq  ,  is  a  function  of 
air  properties  and  orifice  plate  geometry,  static  pressure,  P^  ,  pressure 
drop,  aP  ,  and  temperature  .  Orifice  plate  geometry  was  measured  by  a 
vernier  caliper  having  an  0.4%  error  band.  Static  pressure  was  measured 
from  a  pressure  gauge  with  an  accuracy  of  ±  0.1  psi.  Assuming  an  average 
static  pressure  of  18  psi,  the  resulting  error  is  1.1%.  The  pressure  drop 
was  measured  using  a  water-filled  U-tube  manometer  accurate  to  ±  0.1  in. 
of  water.  Assuming  an  average  AP  of  8  in.  of  water,  this  results  in  an 
error  of  2.5%.  Temperature  was  measured  using  a  standard  thermometer 
accurate  to  1°C.  Assuming  an  average  air  temperature  of  27°C  the  corres¬ 
ponding  error  is  3.7%.  Utilizing  these  individual  error  values  the  net 
error  associated  with  the  measurement  of  air  velocity  is  4.6%. 

Lastly,  the  fluidized  bed  pressure  drops  were  measured  using  an 
inclined  oil-filled  manometer  accurate  to  ±  .01  in.  of  water.  Assuming  an 
average  bed  pressure  drop  of  1  in.  of  water,  there  results  a  2%  error  in 
measurement.  Table  5.1  summarizes  errors  associated  with  all  the  measured 
quantities  discussed  in  this  section. 
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Velocity 

Orifice  Plate 
Static  Pressure 
Pressure  Drop 
Temperature 
Fluidized  Bed 
Pressure  Drop 


TABLE  5.1 

EXPERIMENTAL  ERROR  SOURCES 
ACCURACY 


Wall  Heater 


Tube  Hfcutar 


Heat  Transfer  Rate 
Voltage 
Amperage 
Leakage 
Heater  Area 
Temperature 

Thermocouples 

Difference 


±  1  °F 


±1°F 


I 
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CHAPTER  6 


ANALYSIS  AND  DISCUSSION  OF  EXPERIMENTAL  RESULTS 


This  chapter  will  discuss  in  detail  the  experimental  results  presented 
in  Chapter  5.  Of  concern  in  this  investigation  is  both  the  fluidization 
and  heat  transfer  characteristics  of  the  test  particles  as  compared  with 
the  theoretical  predictions  outlined  in  Chapter  2.  Although  closely 
related,  the  heat  transfer  and  fluidization  subjects  will  be  treated 
separately  in  this  chapter  for  the  purpose  of  clarity.  All  analyses  are 
based  upon  the  detail  presented  in  Chapter  3. 

6,1  Fluidization  Studies 

Fluidization  studies  have  been  conducted  primarily  through  visual 
observation  and  through  measurements  taken  as  described  in  Chapter  5. 

Topics  of  this  section  present  the  fluidization  results  in  terms  of  the 
Gel  dart  classification,  the  presence  of  material  static  electrification, 
fluidized  bed  pressure  drops,  air  pumping  power  requirements,  and  the 
effect  of  particle  diameter  and  density  on  minimum  fluidization  velocities. 
Table  6.1  lists  these  results. 

6.1.1  Geldart  Classi fication 


Chapter  2  presented  the  general  philosophy  behind  Geldart1 s 
materials  classification  according  to  their  fluidization  characteristics. 
This  classification,  however,  was  based  on  materials  significantly  more 
dense  and  generally  of  smaller  diameter  than  the  material  studied  in  this 
work  (see  Figure  2.1).  It  is,  therefore,  of  interest  to  see  whether  the 
general  trends  predicted  by  Geldart' s  classification  can  be  applied  to 
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TABLE  6.1  TEST  MATERIAL  FLUIDIZATION  CHARACTERISTICS 


0) 

(2) 

(3) 

(4) 

p 

T~ 

P 

FT 

mf 

APp 

P/A 

PARTICLE 

[pcf] 

_  [in] 

[fps] 

[psf] 

[HP/ft2] 

• 

Norton  650 

31.8 

0.026 

0.24 

12.71 

0.029 

Macro  M40X 

27.2 

0.039 

1.04 

13.73 

0.136 

Norton  1 

23.2 

0.039 

0.57 

10.61 

0.057 

(* 

Norton  1 .5 

20.9 

0.059 

0.75 

9.44 

0.067 

Macro  M27X 

15.3 

0.078 

1.59 

8.42 

0.127 

Perlite  250 

4.5 

0.010 

0.011 

1.60 

3.2  x  10'5 

• 

Perlite  500 

2.7 

0.020 

0.016 

1.01 

• 

2.97  x  10"5 

NOTES:  (1) 

5p  * 

mean  particle  diameter 

w 

(2) 

■ 

i(u„f 

2-D  Bed 

U  -  ) 

mf 

3-D  Bed 

(3)  AP  = 

i(APB 

apb 

)/ft.  of 

bed  height 

2-0  Bed 

3-D  Bed 

<■> 

(4) 

P/A  = 

Air  pumping 
area  at  0  ^ 

power  per 

unit  bed 

cross-section 

= 

Umf<APB> 

w 


All  quantities  are  experimentally  measured 


larger,  less  dense  materials  in  order  to  be  a  further  aid  in  fluidization 
engineering  or  fluidized  bed  design.  In  this  section  each  test  particle's 
fluidization  characteristics  will  be  discussed  and  compared  to  Gel  dart's 
classification.  Additionally,  the  non-test  particles  described  in  Chapter 
5  will  be  evaluated.  Table  6.2  is  a  summary  of  these  results. 

The  Norton  650  particles  exhibit  very  good  fluidization 
characteristics.  As  shown  in  Fig.  E.l,  they  have  an  average  experimentally 
determined  minimum  fluidization  velocity,  II  ^  ,  of  0.24  fps  and  bed 
pressure  drop  per  foot  of  height,  Apg  of  12.75  psf.  Bubbling  commences 
in  this  material  slightly  above  the  minimum  fluidization  velocity.  Based 
upon  experimental  observations  and  testing,  Norton  650  is  best  described  by 
Geldart's  Group  B  fluidization  classification.  It  also  obeys  Eq.  (2.4) 
which  Geldart  uses  to  describe  a  Group  B  particle.  In  general,  Geldart's 
classification  criteria  satisfactorily  describes  this  material. 

The  Macro  M40X  particles  also  exhibit  good  fluidization 
characteristics.  Figure  E.2  shows  that  this  material  has  an  average 

of  1.04  fps  and  a  APg  of  13.77  psf.  Again,  bubbling  appears  to  commence 
slightly  above  minimum  fluidization.  As  with  the  Norton  650  material,  the 
Macro  M40X  particles  are  best  described  experimentally  by  and  satisfy  the 
analytical  expression  for  Geldart's  Group  B  classification.  Thus,  it  can 
be  concluded  that  Geldart's  classification  also  satisfactorily  describes 
this  material. 

The  Norton  1  material  fluidizes  well.  Bubbles  commence  slightly  above 


minimum  fluidization.  It  has  an  average  and  ApR  of  0.57  fps  and 


B 
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TABLE  6.2 

GELDART  MATERIAL 

CLASSIFICATION 

COMPARISON 

• 

DESCRIPTION 

PARTICLE 

Ap  [kq/m3] 

d  [ym] 

GELDART 

EXPERIMENT 

l 

Norton  650 

507.8 

650 

B 

B 

Macro  M40X 

434.4 

1000 

B 

B 

• 

Norton  1 

370.4 

1000 

B 

B 

(1 

Norton  1 .5 

332.7 

1500 

B 

B 

Macro  M27X 

244.0 

2000 

B/D 

B 

V 

Q-Cel 1 

175.0 

~  70 

A/C 

C 

% 

Perlite  250 

70.9 

250 

A 

A 

Perlite  500 

42.0 

500 

A 

A 

Styrofoam  Beads 

36.4 

3000 

B/D 

C 

% 


10.64  pcf,  respectively,  as  illustrated  in  Fig.  E.3.  As  with  the  preceding 
materials,  Geldart's  Group  B  classification  satisfactorily  describes  the 
fluidization  characteristics  of  Norton  1. 

The  next  group  of  test  particles,  Norton  1.5,  also  demonstrates  good 
fluidization  characteristics  showing  bubble  commencement  slightly  above 
minimum  fluidization.  This  material  has  an  average  Um^  of  0.75  fps  and  a 
APg  of  9.48  pcf  (Fig.  E.4).  Its  fluidization  characteristics  are  well 
described  by  the  Group  B  description  and  it  satisfies  Eq.  (2.4).  Again, 
the  Geldart  classification  provides  a  satisfactory  description  of  this 
material . 

The  Macro  M27X  material  fluidizes  fairly  well,  although  not  as  well  as 
the  preceding  test  materials.  This  material,  primarily  due  to  its  particle 
diameter,  requires  the  greatest  amount  of  air  velocity  for  fluidization 
(  1.59  fps).  As  a  consequence  of  this,  bubble  diameters  appear  to  be 
larger  and  the  bed  fluidization  is  more  violent.  However,  due  to  its  lower 
particle  density,  this  material  has  a  lower  APg  than  the  preceding  test 
materials  (  8.45  pcf)  as  shown  in  Fig.  E.5.  This  material  is  not  very  well 
described  by  Geldart's  materials  classification.  It  appears  to  be 
somewhere  between  the  Group  B  and  D  categories. 

On  the  other  hand,  both  the  Perlite  250  and  Perl  it  500  materials 
fluidize  well.  They  exhibit  some  bed  expansion  before  the  onset  of 
bubbling.  The  Perlite  250  has  an  average  U  ^  of  0.011  fps  and  a  APg  of 
approximately  1.60  pcf  (Fig.  E.6).  The  Perlite  500  has  an  average  of 

0.016  fps  and  a  APg  of  1.02  psf  (Fig.  E.7).  Their  fluidization  is 
hampered,  though,  by  their  tendency  to  attain  and  hold  a  triboelectric 
charge.  This  is  most  apparent  when  the  particles  are  in  contact  with  the 


plexiglas  front  walls  of  the  fluidized  beds.  As  discussed  in  Appendix  C, 
both  perlite  and  plexiglas  are  dielectrics  and  thus  have  relatively  small 
susceptabi 1 i ties  available  for  this  triboelectric  charging  to  take  place. 
Perlite  in  contact  with  the  metal  and  wooden  walls  does  not  present  this 
problem.  The  attraction  between  the  perlite  and  the  plexiglas  is  not 
enough  to  seriously  hamper  the  perlite's  fluidization  character!' sties,  but 
it  may  pose  a  heat  transfer  problem.  Perlite  250,  moreover,  tends  to 
exhibit  a  greater  attraction  to  the  plexiglas  than  the  Perlite  500. 

Despite  this  triboelectrification  problem,  though,  both  perlites  are 
best  categorized  by  Geldart's  Group  A  description  (but  graphically  appear 
in  Group  C).  Thus,  it  appears  that  his  materials  classification  narrative 
descriptions  can  also  apply  to  this  category  of  material. 

In  summary,  then,  Geldart's  material  classification  according  to 
fluidization  characteristics  seems  to  hold  adequately  for  the  lower  density 
and  particle  diameter  materials  used  in  this  investigation.  His  initial 
criteria  based  on  ranges  of  particle  diameters  and  densities  for  the  four 
fluidization  groups,  however,  needs  to  be  modified  to  include  the 
appropriate  diameters  and  densities  of  these  test  materials. 

Before  concluding  this  topic,  though,  a  discussion  of  two  materials 
initially  investigated  in  this  work  whose  fluidization  character!' sties  are 
not  well  described  by  Geldart's  material  classification  is  in  order.  These 
two  materials  are  the  styrofoam  beads  and  the  Q-cell  [27]  (small  diameter, 
moderate  density)  material  described  in  Chapter  5. 

Considering  the  material  density  and  particle  diameter  of  the 
styrofoam  beads  and  Eqs.  (2.4)  and  (2.5)  one  would  predict,  based  upon 
Geldart's  classification,  that  this  material  would  be  in  the  Group  B/D 
range,  being  very  close  to  Group  D.  In  experimental  investigation. 
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however,  this  material  exhibited  some  bed  expansion,  displayed  excessive 
triboelectric  charging  and  in  general  fluidized  poorly.  This  fluidization 
characteristic  is  best  described  as  a  Gel  dart  Group  C  material  and  not  the 
predicted  Group  0. 

The  second  unpredicted  material  is  the  "Q-Cell"  material  of  mean 
particle  diameter  of  70  pm  and  density  of  11  pcf.  Given  this  data,  one 
would  predict  the  material  to  fall  generally  into  the  Group  A  region  due  to 
the  mean  particle  diameter  criteria.  In  reality,  however,  this  material 
fluidizes  very  poorly  and  seems  to  be  better  represented  by  Group  C.  It  is 
interesting  to  note  that  if  the  classification  of  this  material  was  based 
on  Figure  2.1  alone,  the  Group  C  classification  would  be  accurately 
predicted.  If,  however,  the  classification  was  based  on  particle  density 
and  diameter,  without  reference  to  Figure  2.1,  a  Group  A  classification 
would  be  predicted.  Figure  6.1  shows  these  test  results  along  with  those 
previously  reported  by  Geldart. 

Thus,  it  appears  that  care  must  be  exercised  in  using  the  Geldart 
materials  classification  on  the  very  extremes  of  the  low  density  region, 
that  is,  at  the  very  large  particle  diameter  (  2500  pm)  and  the  very  small 
diameter  (  50  ym)  areas  and  when  fluidizing  materials  where  conditions 
exist  for  the  creation  of  significant  triboelectric  charging. 

6.1.2  Air  Pumping  Power  Requirements 

This  section  is  of  interest  due  to  the  intended  applications  of 
this  work  in  solar  collector  and  heat  exchanger  technology.  The  analysis 
presented  here  will  be  further  utilized  in  Section  6.2  and  in  Chapter  7. 

The  external  air  pumping  power  required  to  operate  a  fluidized  bed  is,  in 


itself,  not  a  fluidization  characteristic  of  a  material  but  rather  a  result 


of  certain  characteristics. 
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In  general,  air  pumping  power,  P_  ,  is  a  function  of  air  mass  flow 
rate,  m  ,  pressure  drop,  AP  ,  and  density,  p 


m  AP 


(6.1) 


Defining  m  =  UQA  and  substituting  into  and  rearranging  terms  in  Eq. 
(6.1)  yields 


P/A  =  U0APb 


(6.2) 


which  expresses  air  pumping  power  per  unit  cross-sectional  area. 
Substituting  Eq.  (2.2)  into  Eq.  (6.2)  results  in  an  expression  for  IP/A  in 
terms  of  fluidized  bed  characteristics  at  UQ  >_  11^ 


P/A  =  U0L(1  -  c)(ps  -  cg)  f- 


(6.3) 


Consequently,  either  Eqs.  (6.2)  or  (6.3)  may  be  used  to  express  the  air 
pumping  power  requirement  per  unit  bed  cross-sectional  area. 

As  expressed  by  Eq.  (6.2),  (P/A)  is  linearly  proportional  to  the 
product  of  the  air  velocity  through  the  bed  and  the  bed  pressure  drop. 
Thus,  in  order  to  minimize  (P/A)  for  solar  collector  and  heat  exchanger 
applications,  this  product  must  be  minimized.  Air  velocity  may  be 
minimized  to  slightly  above  minimum  fluidization  velocity.  Therefore, 
materials  with  a  small  U  ^  will  result  in  a  lower  (P/A)  .  The  material 
parameters  affecting  are  discussed  in  Section  6.1.3. 

Bed  pressure  drop,  as  shown  in  Eq.  (2.2),  is  strongly  affected  by  the 
material  density.  Generally  speaking,  lower  density  material  will  yield  a 


smaller  APg  .  Thus,  low  density  material  will  also  aid  in  yielding  low 
(£/ A)  values.  Table  6.1  shows  this  trend  for  (P/A)  values  evaluated  at 
Un  =  Um-  . 

U  m  i 

Consequently,  desirable  materials  for  low  (P/A)  values  should  have 
both  low  U  ^  and  density  characteristics.  As  shown  in  Table  6.1,  the 
Norton  650  and  the  Perlites  in  this  work  require  the  lowest  (P/A)  value. 

6.1.3  The  Effect  of  Particle  Diameter  and  Density  on  Um^ 

In  Chapter  2  it  was  shown  that  Kunii  and  Levenspiel,  using  the 
correlation  of  Eq.  (2.3),  have  analytically  expressed  in  terms  of, 

primarily,  particle  diameter  and  density.  They  stress  that  this  relation 
is  only  approximate  and  valid  for  uniformly  sized  material.  The  test 
particles  in  this  work  are  not  of  constant  particle  diameter  and  this 
equation  cannot  be  applied  directly.  However,  this  expression  is  useful  in 
interpreting  data. 

Specifically,  Eq.  (2.3)  shows  that  is  a  function  of  (dp)  and 

o$  .  Thus,  solely  based  on  this  relation,  one  would  expect  particle 
diameter  to  have  a  stronger  effect  on  than  particle  density.  As 

shown  in  Table  6.1,  this  is  generally  the  case  with  the  test  materials  of 
this  investigation.  The  Macro  M27X  material  (dp  =  0.078  in, 
p  =  15.3  pcf)  has  a  Um^  of  1.59  fps  versus  the  Macro  M40X  (dp  =  0.039  in, 
p=  27.2  pcf)  with  a  of  1.04  fps  demonstrates  this  effect.  However, 

again  care  must  be  taken,  for  this  trend  is  not  an  absolute  as  illustrated 
by  comparing  the  Macro  M40X  and  the  Norton  1.5.  It  appears  that  at 
relatively  close  particle  diameters  density  becomes  significant  in 
determining  Umf  . 

Using  this  conclusion,  lower  U  ^  values  should  correspond  to  smaller 
particle  diameters.  Since,  as  discussed  in  Section  6.1.2,  the  (P/A) 


values  required  to  fluidize  a  bed  of  particles  is  of  concern  in  this  work, 
a  minimization  of  this  value  will  also  result  in  using  smaller  diameter 
materials.  Consequently,  the  best  possible  combination  of  properties  that 
a  material  could  have  in  order  to  minimize  air  pumping  power  is  a  low 
density  and  small  particle  diameter.  In  actuality,  though,  this  criteria 
presents  a  problem.  In  the  manufacturing  of  low  density  materials, 
particle  diameter  is  usually  varied  by  expanding  the  material  with  air  or 
another  gas.  As  a  result,  for  a  specific  group  of  materials  lower  density 
corresponds  to  larger  particle  diameter.  This  is  identically  the  case  for 
all  the  test  particles  in  this  work  for  small  (P/A)  .  This  fact  also 
makes  simple  material  selection  based  on  low  diameter  and  density 
difficult.  It  is  necessary  to  test  various  combinations  of  diameter  and 
density  to  determine  the  minimum  value.  In  this  work,  for  example.  Perlite 
500  has  the  lowest  (P/A)  value  (at  UQ  =  Um^)  even  though  it  does  not 
have  the  smallest  particle  diameter.  Table  6.3  compares  the  theoretical  to 
experimental  values  of  the  test  materials  Um^  and  APg  . 

6.2  Heat  Transfer  Studies 

Heat  Transfer  studies  of  the  test  material  have  been  conducted  for  the 
wall  flat  plate,  immersed  tube,  and  finned  immersed  tube  cases  in 
accordance  with  the  procedures  outlined  in  Chapter  5.  The  data  collected 
is  presented  in  Figs.  6.2  through  6.4  and  is  correlated  and  analyzed  in 
this  section  based  upon  the  theory  developed  in  Chapters  2  and  3. 
Additionally,  this  section  will  discuss  the  experimental  results  of 
particle  diameter,  density,  air  velocity  and  bed  obstructions  on  heat 
transfer  coefficients,  and  conclude  with  a  discussion  on  how  heat  transfer 
compares  with  air  pumping  power  for  the  low  density  materials  used  in  this 
investigation. 
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FIGURE  6.2  Test  Material  Wall  Heat  Transfer  Coefficients 
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The  analysis  presented  in  this  section  considers  the  actual 
experimental  heat  transfer  accomplished  for  each  low  density  test  material 
and  compares  those  results  to  theoretical  predictions.  From  this 
comparison  general  trends  and  conclusions  are  developed.  Each  heating 
geometry  is  discussed  separately. 

Wall  Heat  Transfer 

In  Chapter  3  the  specific  theoretical  relationships  applicable  to  this 
work  have  been  outlined.  However,  no  reference  to  the  particle- to-surface 
residence  time,  t  of  Eq.  (2.6)  is  made.  Quantifying  this  value  has 
traditionally  been  difficult.  In  this  work  no  detailed  procedure  has  been 
used  to  determine  the  residence  times.  Its  value  for  use  in  Eq.  (2.6)  for 
this  analysis  is  assumed,  based  on  visual  observation  and  typical  values 
determined  by  other  investigators  [17].  With  this  in  mind,  for  the  wall 
heat  transfer  case,  the  following  particle- to-surface  residence  times  have 
been  used: 


Particle  Group  U/Umf  V  (sec) 

Macro  and  Norton  1.0  -  1.5  0.8 

1.51  -  2.0  0.7 

2.01  -  2.5  0.6 

2.51  +  0.5 

Perl i te  1 .0  -  1 .5  0.5 

1.51  -  2.0  0.4 

2.01  +  0.3 


The  difference  in  times  between  the  two  groups  is  based  primarily  on 
experimental  observation.  It  was  visually  detected  that  the  Perlite 
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materials  exhibited  more  rapid  mixing  during  fluidization  than  did  the 
Macro  and  Norton  materials.  As  a  result,  there  appeared  to  be  a 
significant  decrease  (on  the  order  of  one  half)  in  wall  residence  time  with 
the  Perlite  as  compared  with  the  other  test  material.  Although  the  actual 
reason  for  this  is  not  exactly  known,  it  is  probably  a  result  of  the 
different  fluidization  characteristics  of  Perlite  (Gel dart  Group  A)  versus 
the  other  materials  (Geldart  Group  B).  Although  this  one  observation 
cannot  be  considered  conclusive,  further  and  more  detailed  investigations 
into  the  role  of  material  fluidization  classification  on  particle-to-wall 
surface  residence  time  may  be  appropriate. 

Having  assigned  these  residence  times  to  the  appropriate  ranges  of 
U/Um<:  ,  a  theoretical  wall  heat  transfer  coefficient  has  been  numerically 
determined  for  each  test  material.  These  results  are  listed  in  Table  6.4 
and  graphically  depicted  in  Fig.  6.5.  Also  listed  in  the  table  are  the 
corresponding  experimentally  determined  wall  heat  transfer  coefficients 
from  the  data  presented  in  Figs.  6.2  through  6.4.  These  experimental  wall 
heat  transfer  coefficients  are  with  obstructions  placed  in  the  bed  (see 
Section  6.2.4) . 

As  illustrated  in  Table  6.4  and  Fig.  6.5,  there  exists  a  fairly  good 

correlation  between  the  theoretically  predicted  wall  heat  transfer 

coefficients,  htheor  »  those  tliat  ^ave  been  measured  experimentally, 

h  „  .  However,  some  departures  do  exist, 
exp 

These  departures  from  theory  are  all  in  the  lower  U/U  ^  areas  of  the 
Group  B  materials.  An  actual  comparison  of  only  U/Um^  values  for  these 
departures  reveals  no  specific  limit  or  trend  for  which  all  particles 
exhibit  this  poor  correlation.  However,  there  does  seem  to  exist  a  general 
trend  when  comparing  the  value  of  the  air  velocity  through  the  bubbles 


TABLE  6.4  WALL  HEAT  TRANSFER  DATA 
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TABLE  6.4  (Cont'd.) 


FIGURE  6.5  Test  Material  Wall  Heat  Transfer  Coefficients  -  Experimental  vs.  Theoretical 
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formed  in  each  material,  (IL  -  u  .  This  data  indicates  that  when  the 
(Uq  -  U  values  are  less  than  0.4  fps  all  of  the  Group  B  materials 
exhibit  poor  heat  transfer  correlations  with  theory.  Conversely .when  the 
value  is  above  0.7  fps  all  Group  B  materials  exhibit  good  correlation  with 
theory.  Between  these  two  values  of  (UQ  -  U^)  theoretical  correlation  is 
apparently  dependent  upon  particle  diameter  and  density.  Density  is 
important  when  particle  diameters  are  approximately  equal,  in  which  case 
the  greater  density  material  requires  the  greater  value  of  (Uq  -  Um^)  for 
theoretical  correlation. 

Poor  theoretical  correlation  with  the  actual  results  at  the  low 

(Uq  -  Umf)  values  of  these  low  density,  large  diameter  Group  B  materials 

suggests  that  the  heat  transfer  model  as  expressed  by  Eq.  (2.7)  may  not  be 

accurate  for  this  region.  Air  velocity  affects  the  bed  heat  transfer, 

according  to  the  model,  most  significantly  in  the  emulsion  and  wall  heat 

transfer.  The  emulsion  heat  transfer  coefficient,  h  ,  sees  this  effect 

e 

via  the  particle  residence  time.  As  discussed  earlier,  for  a  given 
fluidized  material  greater  air  velocity  will  yield  greater  particle  mixing, 
and,  thus,  lower  particle- to-surface  residence  times.  If  the  model  is 
assumed  to  be  correct,  then  the  assumed  residence  time  values  used  to 
determine  h0  must  be  inaccurate.  However,  the  residence  times  required 
for  gc. d  correlation  at  the  low  velocity  values  are  generally  longer  than 
what  was  observed  during  experimentation.  These  required  residence  times, 
treq  ,  are  listed  in  Table  6.4. 

Air  velocity  is  also  a  factor  in  the  convective  component  of  hw  . 
However,  this  is  of  secondary  importance  as  c  -npare'i  to  its  effect  on 

due  to  the  characteristically  low  air  veloc’t’--,  -v.lved  with  fluidizing 
this  material.  Also  at  these  air  velocities  ..nvective  component  of 


c 
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is  much  smaller  than  the  conduction  component  which,  in  turn,  is  not 
dependent  upon  air  velocity  at  all. 

Consequently,  it  must  be  concluded  that  at  low  air  velocities  some 
other  physical  mechanism  may  be  taking  place,  other  than  the  packet  renewal 
theory  used  in  this  model,  in  order  to  account  for  the  unpredicted  low 
actual  wall  heat  transfer  coefficients.  One  possibility  of  this  is  that  at 
these  low  air  velocities  particle  packets  may  have  moved  along  the  heated 
surface  for  a  finite  length  of  time  rather  than  simply,  and  rapidly,  making 
contact  and  moving  away.  Experimental  observations  were  made  to  this 
effect,  but  these  extended  residence  times  were  on  the  order  of  1-2 
seconds  for  the  effected  particles  at  these  velocities.  This  could  account 
for  some  of  the  experimental  heat  transfer  coefficent  departures  from 
theory  but  does  not  explain  the  results  of  the  lowest  velocity  Norton  650 
and  Macro  M40X  data. 

As  mentioned  earlier,  the  theory  generally  does  well  in  predicting  the 
wall  heat  transfer  coefficients  except  for  the  departures  discussed  above. 
The  data  indicates  that  the  use  of  the  conduction  Nusselt  Number  of  6  in  h 
fits  the  data  better  than  the  value  of  12  .  Furthermore,  it  shows  that  the 
controlling  resistance  to  heat  transfer  is  primarily  the  emulsion  component 
of  the  surface  resistance.  Again  the  importance  of  air  velocity  and 
particle-to-surface  residence  time  becomes  evident.  Being  the  controlling 
parameter,  hQ  must  be  maximized  in  order  to  maximize  the  overall  heat 
transfer  from  the  wall  to  the  fluidized  material.  To  do  this,  air  velocity 
can  be  increased  or,  in  order  to  not  simultaneously  increase  air  pumping 
power  requirements,  obstructions  can  be  placed  on  the  wall  or  in  the  bed  to 
enhance  particle  mixing  and  reduce  particle  residence  time.  This  concept 
is  discussed  later  in  the  chapter. 


Of  particular  interest  to  the  applications  of  this  work  due  to  its  low 
density,  moreover,  is  the  Perlite  material.  Unlike  the  Group  B  material, 
the  Perlite's  (Group  A)  wall  heat  transfer  coefficients  through  all  ranges 
of  (U/U^)  are  well  predicted  by  the  theory.  Again  the  controlling  heat 
transfer  parameter  is  hg  but  it  plays  an  even  more  dominant  role  in  the 
Perlite  material  than  in  the  Group  B  materials.  This  suggests  the 
increased  importance  in  decreasing  particle  residence  time  to  enhance  heat 
transfer  for  the  Perlite.  Additionally,  heat  transfer  enhancement  may  also 
be  accomplished  by  taking  measures  to  reduce  the  triboelectric  affect 
exhibited  by  the  Perlite  and  the  plexiglas  front  wall.  In  effect,  reducing 
the  static  attraction  should  enhance  the  fluidization  next  to  the  wall, 
and,  thus,  reduce  the  particle-to-surface  residence  time.  It  should  be 
re-emphasized  that  the  triboelectric  phenomenon  present  in  this  situation 
was  present  with  the  Perlite  against  the  plexiglas  only,  not  with  any  other 
surface,  to  include  the  copper  wall  heater  surface. 

In  summary,  then,  the  data  available  for  the  the  low  density,  large 
particles  used  in  this  investigation  indicates  that  the  heat  transfer  model 
represented  by  Eq.  (2.7)  generally  does  well  in  accurately  predicting  wall 
heat  transfer  coefficients.  The  controlling  heat  transfer  coefficient 
tends  to  hg  .  Problems  do  exist,  however,  with  the  Group  B  material  at 
very  low  air  velocities.  This,  in  turn,  suggests  some  other  physical 
mechanism  may  be  at  work  effecting  their  heat  transfer. 

Unfinned  Tube  Heat  Transfer 

As  with  the  analysis  of  the  wall  heat  transfer,  the  unfinned  tube  heat 
transfer  analysis  is  based  upon  the  relationships  outlined  in  Chapter  3. 

The  particle-to-tube  residence  times  are  assumed,  based  on  typical  values 
determined  by  other  investigators  [17],  as  follows: 
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Particle  Group 

”4nf 

tr(s< 

All 

1 .05 

0.6 

1.05  -  1.10 

0.5 

1.11  -  1.99 

0.4 

2.00 

0.3 

Notice  that  in  the  immersed  tube  case  the  residence  times  versus  U/Umf 

apply  to  all  particle  groups.  Contrary  to  the  wall  case,  there  was  no 
experimental  visual  observation  due  to  the  immersed  nature  of  the  tube 
surface  to  suggest  any  changes  in  residence  times  between  the  different 
material  fluidization  groups  was  required. 

Using  these  assumed  residence  time  values,  theoretical  immersed  tube 
heat  transfer  coefficients  have  been  numerically  determined  for  each 
test  material.  The  results  are  listed  in  Table  6.5  along  with  the 
corresponding  experimentally  determined  immersed  tube  heat  transfer 
coefficients.  The  ratios  of  the  individual  experimental  to  theoretical 
heat  transfer  coefficients  versus  U/Umf  are  graphically  illustrated  in 
Fig.  6.6. 

As  shown  in  Table  6.5  and  Fig.  6.6  there  is  a  very  good  correlation 
between  the  theoretical  heat  transfer  coefficient  predictions  and  the 
actual  results  for  all  materials.  Except  for  the  largest  diameter 
particles  (Macro  M27X  and  Norton  1.5)  the  best  data  fit  of  the  theory  is 
with  using  the  conduction  Nusselt  Number  of  6.  For  the  Macro  M27X  and  the 
highest  air  velocity  Norton  1.5  the  conduction  Nusselt  Number  of  12 
provides  a  better  estimate.  This  is  not  to  say,  however,  that  these 
materials  exhibit  greater  heat  conduction.  The  effect  present  here  is 
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probably  due  to  the  increased  role  of  convective  heat  transfer  as  a  result 


of  the  increased  air  velocity  through  these  particles  relative  to  the  other 
particles  tested  brought  on  by  their  larger  size. 

The  data  also  indicates  that  for  the  Macro  and  Norton  materials  there 
is  no  definitively  controlling  resistance  to  heat  transfer  in  the  immersed 
tube  case.  Both  hg  and  h  are  generally  of  the  same  order.  This 
suggests  that  enhanced  heat  transfer  is  best  accomplished  by  increased  air 
velocity  and  the  use  of  smaller  diameter  particles,  not  necessarily  through 
the  use  of  flow  abstructions. 

On  the  other  hand,  the  Perlite  material  data  indicates  that  h&  is 
again  the  controlling  parameter  in  determining  the  overall  heat  transfer 
coefficient.  Thus,  the  same  options  discussed  for  enhancing  wall  heat 
transfer  should  apply  to  the  immersed  tube  case  for  Perlite. 

In  summary,  then,  the  data  available  for  the  low  density,  large 
particles  used  in  this  investigation  indicates  that  the  heat  transfer  model 
represented  by  Eq.  (2.7)  does  well  in  predicting  unfinned  immersed  tube 
heat  transfer  coefficients.  For  the  Norton  and  Macro  materials  there 
appears  to  be  no  dominate  resistance  to  heat  transfer.  The  Perlite 
materials,  however,  have  their  heat  transfer  strongly  controlled  by  hg 
The  assumed  particle  residence  times  also  appear  to  be  valid. 

Finned  Tube  Heat  Transfer 

The  analysis  of  the  immersed  finned  tube  heat  transfer  proceeds 
similarly  to  that  of  the  unfinned  tube  case.  Of  concern  in  these  tests  is 
the  effect  that  the  placement  of  fins  on  an  immersed  tube  will  have  on  its 
heat  transfer  compared  to  the  unfinned  case.  Test  air  velocities  and 
assumed  particle-to-tube  residence  times  for  the  finned  tube  experiments 
have  been  kept  identical  to  the  corresponding  unfinned  experiments. 
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Theoretical  immersed  finned  tube  heat  transfer  coefficients  have  been 

calculated  for  each  test  material  using  the  surface  efficiency  modification 

procedure  discussed  in  Chapter  3.  Table  6.6  presents  these  results  in 

conjunction  with  the  corresponding  experimentally  determined  finned  tube 

heat  transfer  coefficients  and  the  ratios  of  „  to  h.,  . 

exp  theor 

The  results  of  these  experiments  are  not  surprising.  Using  the  same 
assumed  parti  cl e-to-tube  residence  times  as  the  unfinned  case  yields,  in 
general,  very  poor  theoretical  correlation  to  experimental  data. 

Acceptable  correlation  begins  to  exist  (shown  with  a  "*"  in  Table  6.6 )  at 
the  smaller  diameter  particle  level  and  at  the  higher  air  velocity  regions 
of  the  intermediate  diameter  particles.  The  poor  correlation  is  due  to  the 
decreased  particle  mobility  between  the  fins  in  the  larger  particle  cases. 
The  decrased  particle  mobility  translates  directly  to  an  increased 
parti  cl e-to-tube  residence  time  than  assumed  and,  thus,  a  lower  value  of 
hg  than  would  be  predicted.  The  Perlites  show  good  theoretical 
correlation  because  their  average  diameter,  dp  ,  is  so  small  compared  to 
the  fin  spacing  X  ,  that  the  presence  of  the  fins  in  this  case  is 
unnoticeable.  For  Perlite  500  this  ratio  is  0.06  whereas  for  Norton  1.5, 
it  is  0.18.  Consequently,  the  residence  times  between  the  finned  and 
unfinned  tubes  remain  similar.  Figure  6.7  illustrates  the  effect  of 
increased  particle  diameter  on  the  average  heat  transfer  coefficients  at 
the  given  fin  spacing  of  this  investigation. 

Although  the  model  of  heat  transfer  represented  by  Eq.  (2.7)  is 
difficult  to  apply  to  the  finned  tube  case  due  to  the  uncertainty 
surrounding  the  particle  residence  times,  another  aspect  of  the  theory 
discussed  in  Chapter  2  does  hold  well  in  this  investigation. 
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This  concept,  suggested  by  Glicksman,  suggests  reducing  the  unfinned 


tube  heat  transfer  coefficient  by  30-50i  to  estimate  the  finned  tube 
coefficient.  This  "rule  of  thumb"  works  well  for  these  test  particles. 
Violations  of  this  rule  appear  to  occur  only  where  the  effect  of  the  fins 
on  fluidization  is  small.  In  general,  the  experimental  data  for  these 
materials  indicate  that  where  theoretical  correlation  is  not  good  this  rule 
of  thumb  holds. 

As  would  be  expected,  then,  as  particle  diameter  lengths  approach  that 
of  the  fin  spacing  of  a  finned  surface,  heat  transfer  coefficient  values 
decrease  in  a  fluidized  bed.  However,  the  intended  purpose  of  using  a 
finned  surface  over  an  unfinned  one  is  to  enhance  the  overall  rate  of  heat 
transferred  to  a  medium.  For  the  particular  materials  and  finned  tube 
heater  used  in  this  investigation,  this  enhancement  was  significant.  The  • 
degree  of  enhancement  can  best  be  illustrated  by  comparing  the  unfinned 
tube  heat  transfer  coefficients  to  a  psuedo-finned  tube  heat  transfer 
coefficient  defined  as: 


h' 


«f 

A  ATf 


(6.4) 


where 

Qf  =  rate  of  heat  transfer  for  the  finned  tube  experiment  [Btu/hr] 
A  =  the  tube  surface  area  of  the  finned  tube  heater  without  fins 
[0.074  ft2] 

AT ^  =  the  temperature  difference  between  the  bare  surface  of  the 
finned  tube  heater  and  the  fluidized  material. 


-94- 


By  comparing  hUN  to  h'  gives  an  indication  of  the  heat  transfer 
enhancement  due  to  the  presence  of  fins  is  available.  Table  6.7  shows  this 
comparison.  As  is  indicated  in  this  table,  heat  transfer  has  been  enhanced 
by  a  factor  of  4  to  9  while  the  actual  unfinned  tube  heat  transfer 
coefficient  has  only  been  reduced  by  1/3  -  1/2  for  the  larger  particles  and 
no  more  than  by  1/5  for  the  smaller  particles.  These  results  suggest  that 
an  optimum  particle  diameter  to  fin  spacing  ratio  and  fin  efficiency  may 
exist  in  which  the  reduction  of  hUN  to  hpj N  is  minimized  while 
maximizing  the  heat  transfer  enhancement. 

In  summary,  then,  the  data  available  for  the  low  density,  large 
particles  used  in  this  investigation  indicates  that  the  heat  transfer  model 
represented  by  Eq.  (2.7)  does  well  in  predicting  finned  immersed  tube  heat 
transfer  coefficients  only  in  situations  where  the  ratio  of  particle 
diameter  to  fin  spacing  is  very  small.  As  this  ratio  exceeds  about  0.08 
particle  fluidization  between  the  fins  is  adversely  affected,  residence 
times  increase,  and  heat  transfer  coefficients  decrease.  How  the  particle 
residence  time  is  increased  is  difficult  to  quantify.  Until  it  can  be 
better  evaluated  the  model  cannot  be  accurately  used.  However,  it  has  been 
noticed  that  in  the  cases  where  the  model  cannot  be  applied,  the  reduction 
of  hUN  to  hpIN  generally  follows  the  "30/50  Rule  of  Thumb".  Lastly, 
for  all  test  particles  the  actual  rate  of  heat  transferred  was 
significantly  greater  using  fins  than  the  reduction  of  the  heat  transfer 
coefficient. 

6.2.2  The  Effect  of  Particle  Diameter  and  Density  on  Heat  Transfer 
The  nature  and  variety  of  the  test  material  used  in  this 
investigation  makes  it  difficult  to  set  forth  specific  conclusions  on  the 
effect  changes  in  a  particle's  diameter  and/or  density  would  have  on  heat 
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TABLE  6.7 

FIN  ENHANCEMENT  OF 

HEAT  TRANSFER 

PARTICLE 

hun 

h' 

h'/h 

Norton  650 

15.70 

98.07 

6.2 

16.78 

117.42 

7.0 

17.69 

123.40 

7.0 

18.37 

131.67 

7.2 

Macro  M40X 

11.31 

50.17 

4.4 

12.08 

59.46 

4.9 

13.64 

64.66 

4.7 

13.90 

78.31 

5.6 

Norton  1 

12.90 

70.54 

5.5 

13.28 

77.99 

5.9 

13.55 

91.00 

6.7 

14.47 

107.02 

7.4 

Norton  1.5 

11.22 

50.10 

4.5 

12.50 

60.75 

4.9 

12.75 

76.18 

6.0 

13.86 

87.50 

6.3 

Macro  M27X 

12.20 

56.63 

4.6 

13.03 

62.38 

4.8 

13.73 

72.56 

5.3 

13.37 

73.92 

5.5 

Perlite  250 

7.09 

45.32 

6.4 

7.09 

53.55 

7.6 

7.28 

57.50 

7.9 

7.69 

57.85 

7.5 

Perlite  500 

5.67 

48.09 

8.5 

5.86 

51.16 

8.7 

5.92 

53.28 

9.0 

6.12 

56.00 

9.2 

NOTE:  Units 

of  h  are 

Btu/hr-ft2-°F 
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transfer  coefficients.  To  do  this  would  require  the  ability  to  vary  a 
specific  material's  particle  diameter  while  holding  its  density  constant 
and  vice  versa.  As  was  explained  earlier,  this  is  not  possible  with  the 
materials  on  hand.  Consequently,  only  a  general  view  can  be  extracted  from 
the  results  of  this  investigation.  Figures  6.8  and  6.9  show  the  average 
trends  exhibited  by  the  test  materials  in  a  well  fluidized  bed  on  heat 
transfer  with  varying  particle  diameter  and  density.  In  general  it  can  be 
seen  that  increased  particle  density,  represented  in  hg  ,  appears  to  have 
a  more  significant  effect  on  increased  heat  transfer  than  increased 
particle  diameter,  represented  in  hw  .  That  is,  hgxp  appears  to  change 
strongly  with  hg  ,  but  only  changes  significantly  with  hw  when  the  ratio 

»f  Vhe*p  =  ’-2' 

This  conclusion  should  not  be  surprising.  As  discussed  in  Chapter  2 
and  illustrated  earlier  in  this  chapter,  particle  diameter  has  its 
strongest  effect  on  fluidization  velocities.  Heat  transfer,  on  the  other 
hand,  is  directly  affected  by  the  square-root  of  the  particle  density,  as 
shown  in  Eq.  (2.6).  Thus,  changes  in  density  should  have  a  greater  impact 
on  heat  transfer  than  changes  in  particle  diameter. 

6.2.3  The  Effect  of  Air  Velocity  on  Heat  Transfer 

As  illustrated  in  the  data  previously  presented,  the  amount  of 
air  velocity  through  a  fluidized  bed  is  an  important  factor  for  heat 
transfer.  In  general,  it  would  be  expected  that  as  the  air  velocity  is 
increased,  that  is  as  U/Um^  increases,  the  heat  transfer  coefficient 
would  also  increase.  However,  for  the  materials  used  in  this  investigation 
there  are  various  agreements  and  disagreements  with  this  expectation.  To 
evaluate  this,  a  measure  of  the  air  velocity  effect  on  changes  in  heat 
transfer  coefficients  is  accomplished  by  determining  the  ratio  of  the 
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material  ' s  average  heat  transfer  coefficient,  ,  to  its  maximum  value, 

hMAX  ,  and  comparing  it  to  the  material's  average  U/Umf  value.  This 
nondimensional  comparison  gives  the  relative  magnitude  of  the  effect  that 
changes  in  U/Umf  have  on  changes  in  the  heat  transfer  coefficient.  The 
results  of  this  comparison  for  each  test  material  and  heating  geometry  is 
shown  in  Fig.  6.10. 

Evaluating  the  data  in  this  light  yields  several  conclusions.  First, 

without  exception,  all  test  materials  exhibit  a  more  significant  change  in 

wall  heat  transfer  as  U/Um^  changes  than  for  the  unfinned  tube  heat 

transfer.  In  the  range  of  average  U/Um^  values  used  in  this  work  the 

average  wall  heat  transfer  coefficient  varied  from  approximately  0.65  to 

0.85  of  the  maximum  value.  In  the  unfinned  tube  case,  however,  the  average 

heat  transfer  coefficient  only  varied  from  approximately  0.9  to  0.95  of 

li  „  .  These  results  indicate  the  increased  sensitivity  the  wall  heat 
max 

transfer  coefficient  has  with  air  velocity  over  the  unfinned  tube  one.  Thus 
it  would  appear  that  by  increasing  U/Umf  the  wall  heat  transfer 
coefficient  will,  in  general,  increase  more  than  the  corresponding  unfinned 
tube  coefficient. 

Secondly,  finned  tube  heat  transfer  coefficients  appear  to  exhibit 
behavior  between  that  of  the  wall  coefficients  and  the  unfinned  tube  ones 
as  U/Umf  changes.  Additionally,  in  the  two  immersed  tube  cases, 
increases  in  U/Um^  appear  to  have  a  more  dramatic  effect  on  changes  in  h 
for  materials  of  1  mm  and  greater  diameter  than  in  the  smaller  materials. 

Finally,  the  Perlite  materials  demonstrate,  on  average,  less 
sensitivity  of  h  to  Increases  In  U/Umf  for  all  heating  geometries  than 
do  the  other  materials.  This  is  especially  evident  in  the  immersed  tube 
cases.  Thus,  to  improve  the  heat  transfer  coefficients  for  Perlite 


requires  significant  Increases  in  U/U„r  . 


Increasing  air  velocity  through  a  fluidized  bed  of  these  test 
materials  to  increase  their  heat  transfer  is  not  always  appropriate, 
however,  as  indicated  by  the  data.  All  of  these  materials,  in  both  the 
unfinned  and  finned  tube  experiments,  have  shown  a  tendency  of  either 
approaching  an  asymptotic  value  for  the  heat  transfer  coefficient  or 
reaching  a  peak  value  and  then  beginning  to  decline  as  U/Um^  was  further 
increased  (see  Figs.  6.2  -  6.4). 

Thus,  the  data  for  this  investigation  suggests  several  things. 
Increases  in  U/Um^  generally  yield  greater  increases  in  wall  heat  transfer 
coefficients  than  immersed  tube  coefficients.  Tube  heat  transfer 
coefficients  of  the  larger  diameter  materials  are  more  responsive  to 
increases  in  U/Um^  than  the  smaller  materials.  The  Perlite  material  is 
generally  less  affected  by  increased  U/Umf  than  the  other  materials,  and 
all  materials  have  demonstrated  a  peaking  characteristic  in  the  value  of 
their  immersed  tube  heat  transfer  coefficients. 

6.2.4  The  Effect  of  Bed  Obstructions  on  Wall  Heat  Transfer 

The  previous  topic  discussed  the  effect  increases  in  air  velocity 
have  on  the  heat  transfer  coefficients  of  the  materials  tested  in  this 
investigation.  In  general,  it  has  been  shown  that  increasing  air  velocity 
aids  in  increasing  heat  transfer.  However,  there  is  an  apparent  limit  to 
this  aid.  Excessive  air  velocity  will  lead  to  increased  bed  voidage  and/or 
particle  elutriation,  which  adversely  effects  overall  heat  transfer,  and 
will  lead  to  significant  air  pumping  power  requirements.  One  method 
briefly  investigated  in  this  work  of  enhancing  fluidized  bed  heat  transfer 
without  increasing  air  velocity  is  through  the  use  of  bed  obstructions. 

Obstructions  should  break  up  large  bubble  formations  in  a  fluidized 
bed,  thus  increasing  particle  mixing  and  subsequently  reducing 
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parti  cl  e-to-surf  ace  residence  times.  To  evaluate  this  effect,  one  0.5  inch 
diameter  x  15  inch  wooden  dowel  was  suspended  horizontally  in  the  center  of 
the  fluidized  material  in  the  two-dimensional  bed  used  in  this 
investigation.  Comparison  tests  of  with  and  without  this  obstruction  was 
made  using  the  Perlite  250  and  Norton  1  materials.  The  criteria  used  for 
determining  the  effect  of  the  presence  of  the  obstruction  was  based  solely 
on  resultant  wall  heat  transfer  coefficient  values.  Sample  results  of  this 
comparative  testing  are  shown  in  Table  6.8.  At  low  values  of  U/Umf  the 
presence  of  the  obstruction  had  little  or  no  effect  on  changing  wall  heat 
transfer  coefficients.  However,  at  higher  values  the  heat  transfer 
coefficents  were  improved  by  approximately  25-30%. 

Having  identified  the  value  of  using  flow  obstructions,  it  was 
necessary  to  determine  a  possible  optimum  number  and  positioning  of  the 
obstructions  to  yield  the  greatest  increase  in  wall  heat  transfer.  This 
was  done  using  the  Perlite  500  and  Norton  650  materials.  Yarious 
geometries  of  obstruction  placement  were  tested.  These  geometries  and  the 
test  results  are  listed  in  Table  6.9.  As  a  result  of  these  test  it  was 
concluded  that  obstructions  placed  centrally  to  the  heated  surface  and  in 
the  middle  of  the  bed  performed  better  than  those  placed  either  below  or 
above  the  heated  surface  or  closer  to  either  wall.  For  the  6  inch  high 
wall  heater  used  in  this  work  it  was  further  found  that  two  parallel 
dowels,  one  up  two  inches  from  the  heater  bottom  and  the  other  down  two 
inches  from  the  top,  resulted  in  higher  heat  transfer  coefficients  than  any 
other  dowel  combination. 

The  tentative  conclusions  that  can  be  drawn  from  this  are  that 
obstructions  of  air  flow  through  a  fluidized  bed  appear  to  improve  wall 
heat  transfer  at  moderate  to  high  values  of  U/U  -  . 


Obstructions  need  to 


TABLE  6.8  OBSTRUCTION  TESTING  RESULTS 


MATERIAL 
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MATERIAL 


Perlite  500 


Norton  650 


TABLE  6.9  OBSTRUCTION  GEOMETRY  TEST  RESULTS 


GEOMETRY 


1  each  dowel,  on  plexiglas 
bed  wall,  1  inch  up  from 
air  distributor 


1  each  dowel,  on  plexiglas 
bed  wall,  bottom  of  heater 
unit 


1  each  dowel,  center  of  bed, 
center  of  heater  unit 


1  each  dowel ,  center  of 
bed,  center  of  heater 
unit 


2  each  dowels,  center  of 
bed,  center  of  heater 
unit  4  inches  apart 


3  each  dowels,  center 
of  bed,  center  of  heater 
unit  1  inch  apart 


4  each  dowels,  center  of 
bed,  center  of  heater 
unit  1  inch  apart 


h(Btu/hr-ft  -F) 
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be  symmetrically  located  near  the  heated  surface.  The  number  and  perhaps 
geometry  of  the  obstructions  need  to  be  optimized.  Too  few  obstructions 
does  not  yield  the  desired  result  while  too  many  apparently  interfer  with 
the  fluidization  process.  For  these  reasons  all  the  data  ultimately  taken 
in  this  investigation  and  reported  in  this  chapter  for  the  wall  heat 
transfer  was  taken  using  the  two  0.5  inch  dowels  described  above. 

6.2.5  Heat  Transfer  and  Air  Pumping  Power  Requirements 

A  detailed  discussion  of  air  pumping  power  requirements  for  the 
fluidization  of  the  materials  at  the  conditions  of  this  work  was  conducted 
earlier  in  the  chapter.  This  section  will  compare  these  required  pumping 
power  values  to  the  available  heat  transfer  in  the  same  manner  as  is  done 
for  other  forms  of  heat  exchanging  surfaces. 

In  order  to  effectively  make  this  comparison,  air  pumping  power 
requirements  must  be  expressed  in  terms  of  power  per  unit  wall  surface  area 
to  overcome  friction  and  not  in  terms  of  power  per  unit  bed  cross-sectional 
area.  For  the  two-dimensional  bed  used  in  this  investigation,  the  ratio  of 
the  cross-sectional  area  to  the  surface  area  is  0.13.  Figures  6.11  and 
6.12  graphically  depict  the  relationships  between  each  material's  friction 
power  requirement  per  unit  wall  surface  area,  (P/A)  ,  and  the 
corresponding  wall  heat  transfer  coefficient,  hw  based  upon  the  overall 
bed  pressure  drop  with  the  2  bed  obstructions  described  above  .  This 
result,  modified  appropriately,  will  be  utilized  in  Chapter  7  for  making 
comparisons  with  other  heat  transfer  surfaces. 

6.3  Summary  of  Experimental  Results 

As  a  result  of  this  experimental  investigation,  several  results  and 
conclusions  have  been  developed.  All  of  the  low  density  test  materials 
used  in  the  investigation  have  demonstrated  the  capability  of  being 


fluidized  successfully.  The  Geldart  materials  classification  based  on 
fluidization  characteristics  generally  does  well  in  predicting  the  test 
materials  fluidization  behavior,  except  for  the  lowest  density  materials 
where  electrostatic  effects  are  important.  The  data  indicates  that 
particle  diameter  is  the  controlling  factor  in  a  materials  minimum 
fluidization  velocity  and,  ultimately,  external  air  pumping  power 
requirements.  Density,  on  the  other  hand,  controls  the  fluidized  bed 
pressure  drop. 

Further-more,  the  model  of  heat  transfer  depicted  by  Eq.  (2.7)  does 
reasonably  well  in  predicting  the  overall  heat  transfer  coefficient  from  a 
heated  surface  to  the  fluidized  particles.  At  low  values  of  (UQ  -  Umf), 
however,  the  data  indicates  that  other  mechanisms  may  be  taking  place  than 
those  represented  by  the  model.  Also,  the  model  cannot  be  confidently 
applied  to  a  finned  surface  due  to  the  uncertainty  of  the  effect  the 
presence  of  the  fins  has  on  the  material's  fluidization.  In  this  case, 
reducing  the  unfinned  heat  transfer  coefficient  by  30-50X  for  these 
materials  works  reasonably  well  in  predicting  their  finned  heat  transfer 
coefficients  assuming  the  fin  spacing  is  not  too  small. 

Data  on  the  various  heating  geometries  indicates  that,  in  general,  at 
a  given  fluidization  velocity  the  immersed  tube  heat  transfer  coefficient 
is  larger  than  the  corresponding  one  from  the  fluidized  bed  wall.  However, 
increasing  air  velocity  and  using  flow  obstructions  tends  to  increase  the 
wall  heat  transfer  coefficients  at  a  greater  rate  than  the  immersed  tube 
ones.  The  controlling  component  of  these  overall  heat  transfer 
coefficients  is  the  emulsion  component,  hg  .  Particle  density  is  shown  to 
have  a  significantly  larger  effect  on  heat  transfer  than  particle  diameter. 
Increased  air  velocity  through  the  bed  increases  heat  transfer,  but 


indications  are  that  there  is  an  upper  limit  to  this  before  increased 
voidage  and  particle  elutriation  effects  cause  a  decrease  in  heat  transfer. 

Finally,  and  of  special  importance  to  the  intended  applications  of 
this  work,  the  low  density  test  materials  showed  fairly  good  ratios  of 
energy  delivered  to  external  operating  energy  required.  Chapter  7  will 
discuss  this  and  other  aspects  of  utilizing  these  experimental  results  and 
conclusions  for  specific  engineering  applications. 


CHAPTER  7 


SELECTED  ENGINEERING  APPLICATIONS 


This  chapter  will  apply  the  results  and  identified  trends  of  low 
density  particle  fluidization,  as  discussed  in  Chapter  6,  to  two  specific 
engineering  concepts:  the  use  of  a  two-dimensional  fluidized  bed  as  an 
active  flat  plate  solar  collector  and  the  use  of  a  fluidized  bed  as  a  heat 
exchanger.  Of  primary  concern  in  these  applications,  and  which  will  serve 
as  the  basis  of  this  discussion,  is  the  thermal  performance  and  associated 
costs  in  energy  usage  of  the  fluidized  bed  as  compared  with  conventional 
equipment.  The  primary  application  of  this  work,  a  solar  collector,  will 
be  presented  first. 

7.1  The  Fluidized  Bed  Solar  Collector 

The  concept  of  a  fluidized  bed  of  lightweight  particles  used  within 
the  cavity  of  a  building's  south  facing  wall  as  a  vertical  flat  plate  solar 
collector,  as  discussed  in  Chapter  1,  could  not  be  evaluated  uqtil  data  on 
the  fluidization  and  heat  transfer  characteristics  of  the  lightweight 
particles  was  obtained.  As  a  result  of  the  experimental  work  of  this 
investigation  an  initial  evaluation  of  this  concept  can  now  be  conducted. 
From  the  results  of  this  evaluation  further  and  more  defi native  conclusions 
can  be  drawn  on  the  relative  merit  of  the  fluidized  bed  solar  collector 
concept.  In  this  section,  data  drawn  from  the  experimental  work  on  the 
fluidization  and  heat  transfer  of  the  low  density  test  materials  will  be 
utilized  to  compare  the  predicted  performance  and  feasibility  of  a 
"typical"  Fluidized  Bed  Solar  Collector  to  that  of  conventional  liquid  and 
air-cooled  solar  collectors.  To  accomplish  this  a  brief  description  of  the 
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general  characteristics  and  design  of  conventional  flat  plate  solar 
collectors  is  first  given.  Secondly,  an  analysis  of  a  fluidized  bed  as  a 
solar  collector  is  performed.  With  this  information,  comparisons  are  then 
made  between  a  representative  fluidized  bed  solar  collector,  using  the  test 
materials  of  this  work,  and  comparable  conventional  collectors. 

7.1.1  General  Characteristics  and  Design  of  Conventional 


Flat-Plate  Solar  Collectors 


Flat-plate  solar  collectors  in  present  use  require  either  a 
liquid  or  air  as  the  heat  transfer  medium.  A  typical  commercial 
liquid-cooled  collector  panel  has  an  area  of  18  to  32  ft  ,  a  copper 
absorber  plate  with  copper  tubes  manifolded  into  headers,  a  black  chrome 
selective  absorber  surface,  fiberglass  or  foam  insulation,  and  tempered 
glass  covers.  Air-cooled  collectors  differ  from  the  liquid-cooled  ones 
mainly  in  the  design  of  the  fluid  passage  mechanism.  A  continuous  black 
sheet  of  steel  or  aluminum  replaces  the  plate  and  tubes.  A  lower  metal 
sheet  forms  the  bottom  of  the  air  passage.  Inlet  and  outlet  ports  are 
provided  at  the  sides  or  bottom  of  the  unit,  and  metal  screens  are 
optionally  used  as  absorbers  [28]. 

The  thermal  performance  of  any  type  of  solar  collector  can  be 
evaluated  by  an  energy  balance  which  determines  the  amount  of  incoming 
solar  energy,  I£  ,  transferred  to  the  working  fluid  in  the  form  of  useful 
energy,  qu  .  For  a  flat-plate  collector  of  area  Ac  this  energy  balance 


I-A_  T.a$  =  q  +  a,  +  — ^ 
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where  de  /dt  represents  the  rate  of  internal  energy  storage  in  the 

V 

collector  plate  and  is  negligible  for  thin,  metallic  plates. 

For  convenience  in  terms  of  physical  parameters  and  to  allow  direct 
comparison  between  the  thermal  performance  of  collectors  of  varying  designs 
Eq.  (7.1)  may  be  represented  in  the  alternate  form 

%  =  V'C  ‘c  -  uc(Tf  -  Ta)}  (7‘2) 

where  F'  denotes  the  collector  efficiency  factor  and  physically 
represents  the  ratio  of  two  overall  heat  transfer  coefficients  [29].  The 
numerator  is  the  thermal  resistance  between  the  collector  surface  and  the 
ambient  air.  The  denominator  is  the  thermal  resistance  between  the  working 
fluid  and  its  environment.  It  has  a  strong  dependence  on  the  overall 
collector  heat  loss  conductance,  U  ,  and  the  average  heat  transfer 
coefficient  of  the  working  fluid.  It  is  only  slightly  dependent  on 
temperature.  The  collector  efficiency  factor  also  increases  with 
increasing  plate  thickness  and  plate  thermal  conductivity.  Other  new  terms 
in  Eq.  (7.2)  represent  the  temperature  of  the  working  fluid  and  ambient 
air,  and  Ta  .  Additionally,  the  overall  heat  loss  conductance,  U, 
is  primarily  a  function  of  the  number  and  type  of  covers  used,  the  radiant 
emissivity  of  the  collector  surface,  the  effectiveness  of  the  back 
insulation,  and  the  environmental  wind  speed.  Equation  (7.2)  is  limited, 
however,  for  it  only  yields  the  rate  of  heat  transfer  to  the  working  fluid 
at  a  given  point  along  the  collector  plate. 

In  an  actual  collector  the  fluid  temperature  increases  in  the 
direction  of  flow  as  heat  is  transferred  to  it.  Thus,  to  compare  the 
performance  of  a  real  collector  to  the  thermodynamic  optimum,  it  has  been 


convenient  to  define  a  heat-removal  factor,  FR  ,  as  the  ratio  of  the 
actual  rate  of  heat  transfer  to  the  working  fluid  to  the  rate  of  heat 
transfer  at  the  minimum  temperature  difference  between  the  collector  plate 
and  the  environment  [30].  This  thermodynamic  optimum  corresponds  to  the 
condition  of  the  working  fluid  remaining  at  the  inlet  temperature 
throughout  the  collector,  a  condition  which  can  be  approached  at  high 
working  fluid  Reynold's  numbers.  As  a  result,  FR  approaches  F'  ,  as  an 
upper  limit,  with  increasing  fluid  flow  rate.  Thus,  the  rate  of  useful 
heat  transfer,  qu  ,  can  then  be  expressed  as 

%  1  AcFR  1  ^  ‘c  •  Uc<Tf  -  Ta»  (7'3) 


where  Tf  is  now  the  working  fluid  inlet  temperature.  This  form  of  the 
energy  balance  is  frequently  used  for  design  because  the  fluid  inlet 
temperature  to  the  collector  is  usually  known  or  can  be  specified. 

Lastly,  the  collector  efficiency  n  ,  is  defined  as  the  ratio  of  the 
useful  energy  delivered  to  the  total  incoming  solar  energy,  or,  upon  using 
Eq.  (7.3),  as 


n  = 


c  c 


Fr{  ot  -  Uc 
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Ta> 
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For  a  given  design  of  a  flat-plate  solar  collector  FR  ,  at  ,  and  U( 


are  relatively  constant.  Values  of  these  parameters  may  be  determined 
experimentally  and  have  been  estimated  in  the  literature  [29,30], 


Consequently,  as  indicated  by  Eq.  (7.4)  a  straight  line  relationship  exists 


between  n  and  the  variable  (due  to  environmental  conditions) 


<Tf  -  V 


A  fluidized  bed  solar  collector  could  be  designed  similar  to  that 
of  an  air-cooled  flat-plate  collector.  The  difference  being  that  the  air 
would  be  used  to  fluidize  particles  in  a  south  facing  wall  cavity  of  a 
building,  two  to  four  inches  wide,  between  an  external  solar  absorber 
surface  and  an  interior  insulated  wall.  The  wall  cavity  would  be  filled 
with  lightweight  particles  on  the  order  of  those  used  in  this 

investigation.  The  air  would  flow  from  a  distributor  at  the  base  of  the 

cavity,  around  the  particles  and  embedded  heat  transfer  tube(s),  and  be 
collected  in  a  header  at  the  top  of  the  cavity.  Solar  insolation  is 

transmitted  through  the  front  cover  and  absorbed  by  the  collector  surface. 

The  fluidized  particles,  in  turn,  transfer  the  absorbed  energy  from  the 
back  of  the  collector  surface  to  the  embedded  heat  transfer  tube(s).  A 
view  of  a  fluidized  bed  solar  collector  is  shown  in  Figure  7.1. 

As  with  the  more  conventional  liquid-cooled  flat  plate  solar 
collectors,  useful  energy  can  be  removed  from  the  fluidized  bed  solar 
collector  through  a  water  stream  in  the  immersed  heat  transfer  tube(s). 
Additional  energy  could  be  collected  using  the  fluidizing  air  as  a  solar 
energy  collector  medium  similar  to  an  air-cooled  collector  scheme  in  a 
building  space  heating  design.  This  is  limited,  however,  due  to  the  low 
air  velocities  required  for  fluidizing  the  small,  lightweight  particles. 

If  not  used  for  space  heating  purposes,  the  heated  air  could  be 
recirculated  back  through  the  base  of  the  cavity  and  provide  additional 
energy  to  the  immersed  tubes. 

In  order  to  analyze  the  thermal  performance  of  a  fluidized  bed  solar 
collector  Eq.  (7.3)  must  be  further  examined.  The  terms  A  ,  or  I  , 


and  U  ( -  Ta)  are  design  parameters  or  are  determined  from  design 
parameters  and  environmental  operating  conditions  in  the  same  procedure  as 
used  with  a  conventional  flat-plate  collector.  The  evaluation  of  the 
collector  efficiency  factor,  FR  ,  distinguishes  the  performance  of  a 
fluidized  bed  solar  collector  from  that  of  a  conventional  collector.  As 
before,  FR  represents  the  ratio  of  the  actual  rate  of  heat  transfer  to 
the  working  fluid  to  the  rate  of  heat  transfer  at  the  minimum  temperature 
difference  between  the  collector  plate  and  its  environment.  For  a 
fluidized  bed,  expressed  in  terms  of  the  appropriate  resistances  to  heat 
transfer,  this  is  given  by 
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the  thermal  resistance  between  the  collector 
surface  and  the  fluidized  particles 
the  thermal  resistance  within  the  fluidized  bed 
the  thermal  resistance  between  the  bed  particles 
and  the  working  fluid  in  the  immersed  tubes. 


For  absorbing  surfaces  made  of  thin  metal  sheets. 
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where  h. ,  is  the  heat  transfer  coefficient  between  the  back  surface  of  the 


collector  surface  and  the  bed  particles.  In  this  work  it  is  equal  to 


hwall 


The  effective  thermal  resistance  within  a  fluidized  bed  is 
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where  £  is  the  average  path  length  for  heat  flow  from  the  collector 
surface  to  the  immersed  heat  transfer  tube(s)  and  is  the  effective 

thermal  conductivity  of  the  bed.  The  average  path  length,  £  ,  is  a 

function  of  the  bed  design  and  can  be  approximated  by 


£ 
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The  effective  thermal  conductivity,  kgff  ,  is  a  quantity  whose  value  can 
only  be  estimated  at  this  time.  Work  is  presently  being  conducted  to 
accurately  determine  kg^f  for  a  bed  of  small,  lightweight  particles  [17] 
Previous  work  on  beds  of  greater  density  material  has  indicated  that  kg^ 

is  a  function  of  the  amount  of  material  displaced  within  the  bed  due  to 
bubble  motion,  the  effective  density  and  heat  capacity  of  the  particles, 
and  the  height  in  which  particles  are  displaced  within  the  bed  [31] 
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For  heat  transfer  tubes  made  of  thin  metal,  the  resistance  to  heat 


f  1  ow  i  s 


T  £h  A, 


(7.10b) 


where  Eq.  (7.10a)  is  the  resistance  associated  with  an  unfinned  tube  and 
Eq.  (7.10b)  is  the  corresponding  resistance  of  a  finned  tube.  The  term  e 
in  Eq.  (7.10b)  represents  the  overall  surface  efficiency  as  defined  in 
Chapter  4. 

Thus,  as  with  conventional  flat-plate  solar  collectors,  a  fluidized 
bed  collector  efficiency  can  be  represented  by  Eq.  (7.4).  A  straight  line 
relation  would  again  exist  between  n  and  (T*  -T  )/I  for  the  fluidized 
bed  collector. 

7,1.3  Solar  Collector  Comparisons 

Using  the  relationships  developed  in  the  previous  two  sections, 
the  efficiency  of  a  fluidized  bed  solar  collector  can  be  compared  to  that 
of  a  typical  liquid  and  air-cooled  flat  plate  solar  collector  operating 
under  the  same  environmental  and  geometric  conditions.  In  this  section,  a 
standard  liquid  and  air-cooled  flat-plate  solar  collector  will  be  compared 
with  fluidized  bed  solar  collectors  comprised  of  this  work's  low  density 
test  material.  All  collectors  will  bed  designed  to  the  following 
specifications: 

1.  Collector  Surface  Area,  Ac  =  8  ft  x  4  ft  =  32  ft2 

2.  Air-cooled  and  fluidized  bed  collector  cross-sectional 

area,  A =  4  ft  x  4  in  =  1.33  ft2 

3.  81ack  painted,  copper  collector  plate,  e  =  0.95 


4.  Environmental  Wind  Speed  of  8  mph 


5.  Ambient  air  temperature  of  50°F 

6.  No  back  heat  losses 

7.  Vertically  oriented. 

All  collectors  will  be  evaluated  using  one  and  two  glass  covers.  From 
these  specifications  the  following  additional  data  may  be  obtained  [29]: 


r 
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One  Glass  Cover 
Uc  =  1.2  Btu/hr-ft2-  F 
t  $  =  0.91 

a  =  e  * 

ax  =  0.862 


Two  Glass  Covers 
U£  =  0.70  Btu/hr-ft2-  F 
x  $  =  0.83 
0.95 

ax  =  0.786 


Given  these  design  parameters  common  to  all  three  collector  typer.  the 
individual  collectors  may  now  be  designed  and  evaluated.  Table  7.1  lists 
the  additional  details  of  the  liquid-cooled  flat-plate  collector  used  for 
the  comparison. 

An  air-cooled  solar  collector  can  also  be  designed  and  evaluated  using 
standard  techniques  and  the  specifications  listed  above.  Malik  and  Buelow 
[32]  have  surveyed  the  fluid  mechanics  and  heat  transfer  phenomena  in  air 
collectors  and  have  concluded  that  the  heat  transfer  in  a  smooth  air-cooled 
collector  can  be  expressed  as 


St 


0.0192  Re"1/4 
1  +1.22  Re"1/8  (Pr-2) 


(7.11) 


where  St  =  Re  Pr  Nu  ,  the  convective  heat  transfer  coefficient  in  the 
Nusselt  Number  is  based  on  the  unit  collector  area,  and  the  Reynolds  and 


20- 


TABLE  7.1  LIQUID-COOLED  FLAT  PLATE  SOLAR  COLLECTOR  DATA 


Components  :  Copper  tubes  of  0.4  in  ID,  connected  by  a  .02  in 

thick  plate  at  a  center- to-center  distance  of 
6  in.,  7  tubes  total 

:  Working  fluid  -  Water 

:  Working  fluid  flow  rate  -  1.54  lb/sec 

One  Glass  Cover  :  FR  =  0.93 

:  nLCl  =  °-93{*862  "  l*2(Tf  -  Ta)/Ic> 

Two  Glass  Covers  :  FR  =0.98 

:  ri  =  0.98(.786  -  0.7(T,  -  T.  )/I  J 


Nusselt  Numbers  are  based  on  the  unit  hydraulic  diameter.  Using  Eq.  (7.11) 
to  determine  the  convective  heat  transfer  coefficient  of  the  air-cooled 
collector  and  the  relations  developed  in  Section  7.1.1,  the  design  and 
performance  of  the  collector  can  be  completed  and  evaluated.  Table  7.2 
lists  the  additional  details  for  the  air-cooled  flat-plate  collector  used 
in  the  comparison. 

The  fluidized  bed  solar  collector  used  in  the  comparison  will,  in 
addition  to  meeting  the  general  design  specifications,  consist  of  two  4  ft. 
long,  0.75  in.  ID,  finned  immersed  heat  exchange  tubes  symmetrically  spaced 
2.67  ft  apart  within  the  bed.  The  finned  tubes  used  in  this  collector  are 
of  the  same  geometry  as  that  used  in  the  experimental  section  of  this  work. 


Bed  resistances  are  determined  for  each  test  material  based  upon  the 
experimental  results  discussed  in  Chapter  6  and  on  Eqs.  (7.6),  (7.7)  and 
(7.10b).  Table  7.3  lists  the  data  of  each  test  material,  and  an  "optimum" 
material  to  be  discussed  separately,  for  use  in  the  fluidized  bed  solar 
collector  comparison. 

As  can  be  seen  from  this  data,  for  this  particular  design  of  a 
fluidized  bed  solar  collector,  the  tube  resistance  is  the  controlling  heat 
transfer  resistance  when  compared  to  the  wall  resistance.  This  could  be 
changed,  however,  by  possibly  changing  the  geometry  of  the  immersed  finned 
tubes  or  by  adding  more  tubes.  Using  the  model  for  the  bed  resistance  Rg 


expressed  by  Eq.  (7.7)  results  in  this  resistance  being  small  compared  with 
the  other  heat  transfer  resistances  for  the  Norton  and  Macro  materials,  but 
very  large  for  the  Perlites.  This  is  due  to  the  Perlites  low  density  and 
fluidization  velocity. 


c 
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TABLE  7.2  AIR-COOLED  FLAT  PLATE  SOLAR  COLLECTOR  DATA 


Components  :  Air  Flow  Rate  =  0.8  lbm/sec 

i  ReH  =  29760 

:  h  =1.35  Btu/hr-ft2  -  F 

:  Galvanized  Steel  Air  Ducting,  12“  x  12": 
Entrance  -  3  ft.  straight  section, 

90  rectangular  section  elbow, 

1  ft  straight  section  emptying  into  a 
2"  x  48"  plenum  at  the  base  of  the 
air  collector 

Exit  -  Symmetric  with  the  entrace. 

One  Glass  Cover  :  FR  =  0.52 

:  nAa  =  0.52  £.862  -  1  .2(Tf  -  Ta)/Ic} 


< 

Two  Glass  Covers  :  Fn  =  0.65 


:  H'  =  0.65{.786  -  0.7(T*  -  T  )/I  } 

AC2  r  a  c 
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TABLE  7.3  FLUIDIZED  BED  SOLAR  COLLECTOR  DATA 


MATERIAL 

U0  Cfps] 

Rw 

rb 

rt  fr 

1  cover 

FR  2  covers 

Norton  650 

0.6 

.067 

.010 

.151 

0.78 

0.86 

Norton  1 

1.2 

.083 

.008 

.177 

0.76 

0.84 

Norton  1.5 

1.5 

.100 

.007 

.193 

0.74 

0.83 

Macro  M40X 

1.6 

.091 

.007 

.212 

0.73 

0.82 

Macro  M27X 

2.0 

.091 

.018 

.236 

0.71 

0.81 

Perlite  250 

.025 

.105 

1.75 

.321 

0.28 

0.40 

Perlite  500 

.030 

.139 

2.91 

.353 

0.20 

0.30 

"Optimum" 

0.14 

.091 

.157 

.265 

0.53 

0.74 

NOTE:  Units 

for  all  resistances 

are  [hr-ft2-  F/Btu]. 

An  indication  of  the  comparable  performances  of  the  three  types  of 
flat-plate  solar  collectors  can  be  accomplished  by  graphically  depicting 
collector  efficiency,  n  ,  as  a  function  of  (T.  -  T  )/l  .  Using  the 

data  contained  in  Tables  7.1  through  7.3  and  by  applying  Eq.  (7.4),  this 
graphical  comparison  is  shown  in  Figures  7.2  and  7.3.  As  illustrated  in 
these  figures,  the  thermal  performance  of  a  fluidized  bed  solar  collector 
can  be  expected  to  fall  between  that  of  a  comparable  conventional  liquid 
and  air-cooled  collector.  The  exception  to  this  is  in  the  case  of  very  low 
density  fluidized  material,  such  as  the  test  Perlites.  Their  performance 
is  poor  due  to  their  large  bed  resistance  to  heat  transfer.  This  result 
indicates  a  need  to  take  steps  to  reduce  this  resistance  or,  due  to  its 
inherent  uncertainty,  determine  a  more  accurate  measure  of  its  value  to 
allow  for  the  viable  use  of  these  particular  low  density  materials  in  this 
application. 

A  third  option  is  to  use  another  low  density  material,  unavailable  for 
this  investigation,  that  is  of  slightly  greater  density  but  of 
approximately  the  same  size  as  the  Perlite  material.  As  described  in 
Chapter  4,  Perlite's  density  and  particle  diameter  can  be  independently 
controlled  and,  thus,  may  be  a  possible  candidate  for  this  use. 

The  criteria  that  can  be  used  to  determine  the  characteristics  of  this 
"optimum"  material  is  a  direct  result  of  the  analysis  of  Chapter  6.  As 
stated  there,  heat  transfer  showed  a  general  trend  of  increasing  as  density 
and  velocity  increased.  Fluidization  velocity  increased  with  particle 
diameter.  With  these  two  trends,  the  "optimum"  material  should  be  denser 
than  the  present  Perlite  material  and  at  least  as  large  as  the  Perlite  500 
to  yield  greater  heat  transfer.  The  air  pumping  power  required  for 
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FIGURE  7.3  Solar  Collector  Comparisons  (Two  Cover  Glasses) 


fluidization  was  also  strongly  effected  by  material  density  and  the 
required  fluidization  velocity.  Thus,  the  "optimum"  material  should  also 
be  less  dense  and  have  a  smaller  particle  diameter  than  the  Norton  and 
Macro  materials.  The  resulting  suggested  characteristics  of  this  "optimum" 
material  is  a  material  having  a  density  of  approximately  10  pcf  and  a 
particle  diameter  approximately  equal  to  that  of  the  Perlite  500  of  0.02 
in. 

This  material,  perhaps  Perlite-like  but  slightly  denser,  should  have 
the  following  fluidization  and  heat  transfer  characteristics  based  upon  the 
results  of  this  work  and  Eqs.  (2.2),  (2.3),  and  (2.7),  which  have  been 
shown  to  predict  these  characteristics  well: 

Umf  "  4  Umf  Per^lte  ^00  (due  to  density 
increase)  s  0.07  fps 

APg  =  44Pg  of  Perlite  500  (due  to  density) 

=  4  psf/ft  of  bed  height 

hw  a  11  Btu/hr-ft2-°F  at  UQ  =  0.14  fps  (U/Umf  =  2) 

hp  s  8  Btu/hr-ft2-°F  at  Uq  =  0.14  fps  . 

Using  this  data  in  the  same  manner  as  with  the  test  material,  the  fluidized 
bed  solar  collector  data  (Table  7.3)  and  predicted  efficiency  data  (Figures 
7.2  and  7.3)  is  obtained. 

As  indicated  in  Table  7.3  the  increased  density  and  fluidization 
velocity  of  the  ideal  material  has  significantly  lowered  the  bed 
resistance,  Rg  ,  as  compared  with  the  Perlite  500  (the  material  with  equal 
particle  diameter).  Rw  and  Ry  have  remained  on  the  same  order  as  with 
the  test  materials.  Consequently,  the  significant  decrease  of  R„  is 
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solely  responsible  for  the  "optimum"  materials  drastics  improvement  in 
performance  ovei  le  test  Perlite's.  The  additional  importance  of  the 
identification  of  this  "optimum"  material  will  be  presented  later  in  this 
section. 

To  this  point  in  the  discussion  only  thermal  performance  comparisons 
of  the  fluidized  bed  collector  to  a  liquid  and  air-cooled  collector  have 
been  considered.  This  is  of  primary  importance  in  order  to  tentatively 
determine  the  feasibility  of  the  fluidized  bed  solar  collector  concept. 
Having  now  indicated  that  thermally  the  fluidized  bed  solar  collector  has 
the  potential  of  falling  in  the  same  operational  range  as  conventional 
flat-plate  liquid  and  air-cooled  collectors,  other  operational  comparisons 
need  to  be  conducted. 

One  of  these,  the  external  power  requirements,  is  listed  in  Table  7.4 
for  each  solar  collector  based  upon  the  collector  surface  area.  It  is 
informative  to  note  that  the  fluidized  bed  solar  collector  requires  less 
water  pumping  power  than  the  liquid-cooled  collector  for  the  same  water 

p 

heat  transfer  coefficient  inside  the  tubes  (-1000  Btu/hr-ft  -°F).  Also, 

the  power  required  to  pump  the  water  through  a  fluidized  bed  collector  is 

significantly  less  than  that  required  to  fluidize  the  bed, except  for 

Perlite.  All  materials  require  less  air  pumping  power  than  the 

corresponding  air-cooled  collector  in  this  design.  Additionally,  the 

air-cooled  collector  in  this  design,  a  8  ft  x  4  ft  x  4  in  collector  with  a 
2 

1  ft  air  exit  duct  three  feet  long,  is  only  well  suited  for  a  building 
space  heating  scheme  in  which  the  heated  air  is  just  circulated  to  an  open 
space  directly  adjacent  to  the  collector.  The  other  two  collector  types 
deliver  their  usable  heat  in  a  more  versatile  form,  through  a  heated  water 
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TABLE  7.4  SOLAR  COLLECTOR  EXTERNAL  POWER  REQUIREMENTS  [HP/ft*] 


COLLECTOR  TYPE 


WATER 


Liquid-cooled 

2.0  x  10 

Air-cooled 

10  x  10"3 

Fluidized  bed: 

Norton  650 

4.6  x  10‘3 

1.4  x  10 

Norton  1 

7.4  x  10"3 

II 

Norton  1.5 

8.2  x  10"3 

II 

Macro  M40X 

13.0  x  10“3 

II 

Macro  M27X 

9.8  x  10"3 

II 

Perlite  250 

2.3  x  10"5 

II 

Perlite  500 

1.8  x  10"5 

II 

"Optimum" 

3.3  x  10‘4 

II 

stream.  If  the  air-cooled  collector  were  designed  to  deliver  heat  to  the 
building  interior,  its  air  pressure  losses  would  be  even  larger  due  to  the 
more  extensive  ducting  required. 

A  second  operational  consideration  involving  the  fluidized  bed  solar 
collector  is  its  actual  use  in  a  building  wall  cavity.  Generally,  vertical 
walls  are  designed  to  resist  a  laterial  wind  pressure  load,  acting  either 
inward  or  outward,  of  about  20  psf  for  a  wall  less  than  50  feet  high,  plus 
any  additional  loading  the  designer  may  expect  the  wall  to  encounter  [33]. 
If  a  building  wall  cavity  is  to  serve  as  a  fluidized  bed  the  additional 
lateral  load  placed  on  the  wall  due  to  the  bed  pressure  drop  must  be 
considered  in  the  wall  design.  Chapter  6  presented  the  results  of  the  bed 
pressure  drops  associated  with  the  test  materials.  Table  7.5  lists  these 
results  based  upon  an  8  foot  high  wall/fluidized  bed  solar  collector.  As 
shown  in  this  table  the  additional  wall  loading  associated  with  the 
operation  of  a  fluidized  bed  is  significant  for  the  Norton  and  Macro 
materials.  Use  of  these  materials  would,  therefore,  require  a  careful  wall 
design.  The  Perlite  materials,  on  the  other  hand,  yield  little,  if  any, 
additional  wall  loading  over  what  is  generally  designed  to  accommodate  wind 
effects.  Thus,  the  additional  wall  loading  due  to  the  bed  pressure  drops 
may,  in  itself,  be  either  an  economic  or  construction  limitation  on  the 
ultimate  choice  of  material  to  use  in  a  fluidized  bed  solar  collector.  For 
this  reason  the  importance  of  identifying  the  previously  discussed 
"optimum"  type  material  becomes  evident.  In  a  realistic  context,  the  very 
large  pressure  drops  associated  with  the  Norton  and  Macro  materials, 
despite  their  good  thermal  efficiencies,  could  eliminate  them  for  use  in 
this  application.  Additionally,  the  test  Perlite  materials  are  poor 
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thermal  performers.  The  "optimum"  material  provides  a  compromise,  good 
thermal  and  wall  load  results. 

Several  other  considerations  are  also  involved  in  the  comparison  of  a 
fluidized  bed  solar  collector  to  conventional  ones.  The  fluidized  bed 
solar  collector  cannot  be  tilted  at  large  angles  for  the  vertical  as  can 
liquid-cooled  collectors  in  order  to  increase  the  amount  of  solar 
irradiation  received  unless  multiple  air  distributors  are  used. 

Furthermore,  the  wall  cavity  containing  the  bed  must  be  carefully  sealed  to 
prevent  moisture  or  other  foreign  matter  from  entering  the  bed  and  to 
prevent  the  loss  of  bed  particles. 

Conventional  flat-plate  col!ec+ors,  on  the  other  hand,  despite  their 
proven  performance  and  versatility  have  had  significant  historical 
problems.  Liquid  collectors  have  suffered  from  corrosion,  freeze  damage, 
outgassing  of  volatile  materials,  condensation  beneath  and  between  glass 
covers,  rain  water  leakage,  separation  of  parts,  thermal  distortion, 
charring  and  burning  [28].  Air  collectors  have  suffered  from  some  of  these 
same  problems  and  are  inherently  at  a  disadvantage  due  to  the  low  heat 
transfer  coefficient  of  air  and  the  higher  air  pumping  power  requirements. 

Conventional  collectors  a1 so  offer  little  or  no  insulating  value 
during  periods  of  non-operation.  The  fluidized  bed  collector  would  serve 
as  a  good  insulator  when  defluidized.  As  descried  in  Chapter  4,  the 
materials  used  in  this  investigation  are  all  typically  used  in  lightweight 
or  insulation  applications.  In  the  defluidized  state,  the  effective  bed 
thermal  conductivity  is  primarily  a  function  of  the  entrapped  gas  within 
the  bed.  As  such,  a  defluidized  bed  of  the  Norton  or  Macro  materials  would 
have  a  R  value  of  approximately  6  hr-ft  -°F/Btu  while  Perlite  would  have 


L 


* 
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a  slightly  larger  value.  These  values  compare  very  well  with  typical 

o 

commercial  wall  insulation  R  values  of  3-5  hr-ft  -°F/Btu.  Additionally, 
the  good  insulating  property  of  the  defluidized  bed  would  aid  in  preventing 
freeze-ups  during  power  failures. 

In  summary,  then,  the  experimental  results  of  this  work  indicate  that 
the  use  of  a  fluidized  bed  as  a  fiat-plate  solar  collector  appears  to  be  a 
viable  concept.  Its  thermal  performance  generally  falls  in  between  that  of 
typical  liquid  and  air-cooled  conventional  collectors.  Air  pumping  power 
requirements  are  expected  to  be  less  than  that  required  for  an  air 
collector,  and  water  pumping  power  requirements  less  than  what  is  required 
for  a  liquid  collector.  The  pressure  drops  associated  with  the  greater 
density  Norton  and  Macro  test  materials,  though,  present  a  probable 
structural  problem  for  the  wall  containing  the  fluidized  bed.  During 
off-hours  the  defluidized  bed  would  serve  a  dual  role  as  the  wall  cavity 
insulation. 

It  must  be  re-emphasized,  however,  that  the  comparison  between  a 
fluidized  bed  solar  collector  and  a  conventional  liquid-cooled  and 
air-cooled  flat-plate  collector  conducted  in  this  section  represents  a 
comparison  of  only  one  of  several  possible  vertical,  active  flat-plate 
solar  collector  designs.  The  particular  design  used  in  this  section  was 
chosen  as  a  representative  one  that  also  made  the  best  use  of  the  available 
data  on  the  fluidization  and  heat  transfer  of  low  density  materials.  Other 
possible  designs,  such  as  reducing  the  air  gap  from  4  inches  to  a  smaller 
value  and  using  a  roughened  collector  surface  for  the  air-cooled  collector, 
could  be  considered.  However,  for  an  accurate  comparison  the  fluidized  bed 
solar  collector  design  should  similarly  be  changed  to  reflect  the  reduce 


gap  and  roughened  collector  surface  design  changes.  To  date,  data  on  the 
effect  of  varying  gap  width  and  roughened  wall  surfaces  on  the  fluidization 
and  heat  transfer  of  low  density  materials  is  not  available.  In  the  next 
section,  though,  the  heat  transfer  associated  with  air  flow  over  a  smooth 
and  roughened  heat  exchange  surface  is  compared  with  the  results  of  the 
wall  heat  transfer  of  the  experimental  fluidized  bed  work.  This  comparison 
could  give  some  initial  insight  into  the  effect  of  roughening  an  air-cooled 
collector  surface  and  how  it  compares  with  the  fluidized  bed. 

7.2  Heat  Exchanger  ADDli cations 


Another  area  of  potential  interest  in  the  applications  of  low  density 
material  fluidization  and  heat  transfer  is  in  a  heat  exchanger 
configuration.  This  concept  again  involves  the  use  of  a  building  wall 
cavity  as  a  fluidized  bed.  In  this  application,  however,  the  bed  would 
serve  as  an  energy  diffuser  rather  than  an  energy  concentrator,  as  in  the 
fluidized  bed  solar  collector  application.  The  wall  containing  the 
fluidized  bed  could  serve  a  dual  function  for  use  in  either  a  building 
space  heating  or  cooling  scheme.  In  both  cases,  the  fluidized  bed 
characteristic  of  being  isothermal  during  fluidization  would  be  utilized. 

The  space  heating  scheme  would  consist  of  using  the  fluidized  bed  to 
transfer  heat  from  either  an  outside  solar  collector  plate  (in  the  case  of 
an  exterior  wall),  from  hot  water/steam  passing  through  immersed  heat 
transfer  tubes,  or  from  electric  resistance  heaters.  In  this 


configuration,  the  fluidized  bed  would,  due  to  its  i sothermal 1 i ty,  yield  a 
relatively  constant  temperature  wall  surface.  The  wall  temperature  could 
be  maintained  at  a  comfortable  level  by  controlling  the  incoming  heat  flux 
to  the  bed  and/or  the  bed  operation,  and,  due  to  its  surface  area,  still 


provide  a  significant  heat  flux  to  the  space  to  be  heated.  The 

concentrated  incoming  heat  flux  would  have  to  be  controlled  in  such  a 

manner  to  overcome  the  appropriate  thermal  resistances  at  the  heat  source 

and  through  the  bed.  Thus,  in  this  application  the  critical  resistance  to 

heat  flow  which  effects  the  maintenance  of  a  constant  desired  wall 

temperature  is  at  the  wall  itself,  Rw  . 

Table  7.6  lists  the  wall  heat  transfer  coefficients  (h,  =  1/R  )  of 

w  w 

this  investigation's  test  materials  obtained  from  the  experimental  data. 

In  this  table  a  comparison  is  made  between  the  test  materials  in  a 

fluidized  bed  and  air  flowing  over  a  smooth  flat  plate,  in  the  form  of  the 

air  pumping  power  requirement  per  unit  surface  area,  for  the  same  heat 

transfer  coefficient.  As  can  be  seen  by  these  results,  the  fluidized  bed 

generally  requires  less  power  for  the  same  heat  transfer  than  air  alone; 

the  Perlite  and  optimum  materials  require  significantly  less.  If,  for 

example.  Perlite  250  was  fluidized  in  a  wall  cavity,  the  wall  surface 

temperature  maintained  at  80°F,  and  the  fluidization  velocity  adjusted  to 

2 

yield  a  wall  heat  transfer  coefficient  of  9.5  Btu/hr-ft  -°F,  a  heat  flux  of 
2 

760  Btu/hr-ft  would  be  transferred  from  the  wall  into  the  space  to  be 
heated.  Also,  as  with  the  solar  collector  application,  the  materials 
within  the  wall  when  defluidized  would  serve  as  effective  heat  insulation. 

A  building  space-cooling  scheme  would  operate  similar  to  that  of  the 
space-heating  one.  In  this  case,  however,  when  cooling  is  desired,  such  as 
at  night,  heat  transfer  would  take  place  from  the  heated  interior  wall 
panel,  through  the  bed,  and  to  either  an  immersed  heat  transfer  tube  or  to 
the  cooler  exterior  panel.  Again,  to  insulate  from  the  heat  transfer  the 
bed  would  be  defluidized.  Contrary  to  the  solar  collector  application. 


these  applications  best  suit  the  use  of  Perlite-like  materials  due  to  their 
low  fluidization  power  requirements,  bed  pressure  drops,  and  still 
increased  surface  heat  transfer  coefficients  than  air  flowing  over  a  flat 
plate. 

Another  area  of  interest  related  to  the  heat  exchanger  applications, 
and  mentioned  in  the  last  section,  concerns  the  comparison  of  the  thermal 
performance  of  the  fluidized  bed  to  that  of  air  over  a  roughened  surface. 
Han,  et  al .  [34]  have  conducted  extensive  work  on  the  heat  transfer  and 
friction  power  requirements  of  air  through  rib-roughened  parallel  plates  of 
constant  heat  flux.  Table  7.7  compares  the  results  of  the  fluidized  bed 
wall  heat  transfer  to  that  of  Han's  results  for  air  flowing  at  0.8  lbs/sec 
through  a  roughened  duct  of  the  same  geometry  as  the  air-cooled  solar 
collector  described  in  Section  7.1.  As  indicated  in  these  results, 
roughening  the  surface  over  which  air  flows  increases  its  heat  transfer  and 
also  the  associated  friction  power  requirement.  Under  these  conditions, 
the  air  heat  transfer  coefficient  begins  to  approach  that  of  the  Perlite's 
but  the  power  requirement  is  approximately  two  orders  of  magnitude  greater. 
Thus,  it  would  again  appear  that  use  of  a  Perlite  type  material  in  a 
fluidized  bed,  in  which  only  the  wall  resistance  to  heat  flow  was  of 
concern,  would  generally  perform  better  than  just  using  air. 


TABLE  7.7  A 

COMPARISON  OF 

HEAT  TRANSFER  COEFFICIENT  AND 

POWER 

REQUIREMENT  FOR 

A  FLUIDIZED  BED  AND  AIR  FLOW 

THROUGH  RIB- 

ROUGHENED  PARALLEL  PLATES 

MATERIAL 

h  (  ■) 

hwa11 

(P/A)  (2) 

Uo,3) 

Air  (roughened 

5.2 

3.8xl0~4 

8.0 

surface) 

Norton  650 

14.6 

1 .8xl0"3 

0.67 

Norton  1 

14.4 

4.2xl0_3 

1.7 

Norton  1 .5 

13.0 

4.6xl0"3 

2.1 

Macro  M40X 

15.0 

6.4xl0~3 

2.1 

Macro  M27X 

12.7 

4.9xl0~3 

2.9 

Perlite  250 

9.5 

9. 9x1  O'6 

.02 

Perlite  500 

7.1 

6.7xl0‘6 

.03 

NOTE:  Units  are:  (1)  [Btu/hr-ft2-°F] 
(2)  [HP/ft2] 


CHAPTER  8 


CONCLUSIONS  AND  RECOMMENDATIONS 

In  this  work  the  fluidization  and  heat  transfer  characteristics  of  low 
density  particles  in  a  fluidized  bed  have  been  investigated.  Seven  groups 
of  particles  ranging  in  diameter  from  approximately  250  ym  to  2mm  and 
density  of  2.5  to  32  pcf  have  been  used  for  all  studies.  It  has  been  shown 
that  these  materials  can  be  successfully  fluidized  and  that  their 
fluidization  behavior  can  be  fairly  accurately  predicted  using  the  Geldart 
materials  classification.  Exceptions  to  this  exist,  however,  for  materials 
that  exhibit  significant  triboelectric  effects.  Additionally,  theoretical 
expressions  for  determining  minimum  fluidization  velocities  (Eq.  2.3)  and 
bed  pressure  drops  (Eq.  2.2)  do  well  in  predicting  the  actual  values  for 
these  low  density  materials.  It  has  also  been  shown,  though,  that  not  all 
low  density  particles  can  be  fluidized  for  a  variety  of  reasons 
(i nter-particle  forces,  static-electric  attraction). 

The  model  of  heat  transfer  that  has  been  developed  using  significantly 
denser  materials  than  those  of  this  work  (Eq.  2.7)  has  also  been  shown  to 
do  well  in  predicting  the  overall  heat  transfer  coefficient  from  a  heated 
surface  to  a  bed  of  low  density  fluidized  particles.  However,  it  cannot  be 
applied  with  confidence  to  a  finned  surface  without  knowing  the  effect  the 
presence  of  the  fins  have  on  the  material's  fluidization.  Despite  this 
uncertainty,  though,  the  use  of  fins  does  significantly  increase  the  rate 
of  heat  transfer  to  the  fluidized  bed  and  a  good  estimate  of  the  overall 
heat  transfer  coefficient  can  be  made  using  Glicksman's  "30/50  Rule". 

The  results  of  the  heat  transfer  studies  have  also  indicated  that 
particle  density  is  a  larger  controlling  factor  in  heat  transfer  than 
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particle  diameter.  Density  has  a  significant  input  in  the  emulsion 
component,  hg  ,  of  the  overall  heat  transfer  coefficient,  which,  in  turn, 
has  been  shown  to  be  the  controlling  component.  Particle  diameter,  on  the 
other  hand,  plays  a  dominant  role  in  fluidization  velocities,  but  due  to 
the  inherently  low  velocities  required  for  the  fluidization  of  these 
materials  enters  into  the  heat  transfer  as  a  secondary  contributor. 

The  potential  for  additional  heat  transfer  enhancement  has  also  been 
studied  through  the  use  of  air  flow  obstructions  placed  in  the  bed.  It  is 
hoped  that  the  placement  of  obstructions  would  break  up  large  bubble 
formations  and  thereby  reduce  parti cle-to-surface  residence  times.  Initial 
results  show  that  heat  transfer  coefficients  may  be  increased  up  to  302  at 
moderate  to  high  U/Umf  levels  through  the  use  of  obstructions. 

Using  the  results  of  the  low  density  material  fluidization  and  heat 
transfer  studies,  this  work  has  additionally  investigated  the  feasibility 
of  using  a  fluidized  bed  of  low  density  particles  as  part  of  a  building's 
efficient  energy  management  program.  Two  applications  have  been  discussed: 
integrating  a  fluidized  bed  into  a  building  wall  cavity  to  serve  as  a 
"flat-plate"  solar  collector  and  to  serve  as  a  heat  exchange  medium. 

In  the  solar  collector  application  the  thermal  efficiency  of  a 
fluidized  bed  solar  collector  can  be  expected  to  fall  in  between  that  of 
comparable  conventional  liquid  and  air  cooled  flat  plate  solar  collectors. 
Exceptions  to  this  exist  when  very  low  density  material  (Perlite  of  this 
investigation)  is  used  in  the  fluidized  bed.  This  material  has  a  large 
resistance  to  heat  transfer  through  the  fluidized  bed  due  to  its  low 
density  and  fluidization  velocities  (see  Eqs.  (7.7)  and  (7.9)). 

Additionally,  the  fluidized  bed  solar  collector,  using  all  the  test 


materials,  requires  less  water  pumping  power  than  the  liquid-cooled 


collector  for  the  same  heat  transfer  to  the  water  and  significantly  less 
complex  plumbing  requirements.  Similarly,  the  fluidized  bed  solar 
collector,  using  all  test  materials,  requires  less  air  pumping  power  than 
an  air-cooled  collector  operating  at  normal  conditions  (air  velocity  of  8 
fps)  and  provides  the  added  benefit  of  having  a  more  versatile  and  easily 
transportable  energy  stream  (heated  water  in  a  tube  versus  heated  air  in  a 
duct).  The  fluidized  bed  solar  collector  during  non-operation  also  acts  as 
effective  wall  insulation  in  contrast  to  conventional  collectors. 

However,  the  Norton  and  Macro  materials  used  in  this  work  pose 
potential  wall  structural  problems.  Due  to  their  density  they  create 
significant  pressure  drops  that  may  be  impractical  for  use  in  a  building 
wall  cavity.  Thus,  because  of  the  Perlite's  poor  thermal  performance  and 
the  Norton  and  Macro's  excessive  wall  pressures  an  "optimum"  material 
having  a  density  of  approximately  10  pcf  and  a  diameter  of  0.5  mm  (yet  to 
be  identified)  has  been  described  and  analyzed  for  use  in  this  application. 
This  material,  in  theory,  shows  both  good  thermal  and  wall  pressure 
resul  ts. 

In  the  heat  exchanger  application  the  wall  containing  the  fluidized 
bed  could  serve  a  dual  function  for  use  in  either  a  building  space  heating 
or  cooling  program.  In  this  application  the  use  of  a  Perlite-like  material 
would  be  most  appropriate  due  to  its  low  density  (thus  low  wall  pressures), 
low  fluidization  velocities  (thus  low  air  pumping  power  requirements) ,  and 
because  the  critical  resistance  to  heat  transfer  in  this  application  is  at 
the  wall  surface,  hw  . 

Finally,  Chapter  1  noted  that  little  work  has  previously  been  done  on 
the  use  of  low  density  particles  in  a  fluidized  bed  and  the  potential 
applications  of  such  use.  This  investigation  has  attempted  to  answer  some 
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of  the  questions  surrounding  this  topic.  Although  it  appears  that  low 
density  particle  fluidization  generally  follows  the  same  theoretical 
fluidization  and  heat  transfer  processes  as  greater  density  particles  do 
and  that  feasible  applications  of  this  technology  do  exist,  many  questions 
still  remain  unanswered. 

This  work  has  only  identified  nine  groups  of  low  density  particles  for 
testing  in  which  seven  were  extensively  used.  It  is  necessary  that  future 
work  be  conducted  on  other  low  density  materials  to  develop  further 
experimental  data.  Initially,  it  is  recommended  that  the  earlier  described 
"optimum"  material  be  identified  and  tested.  At  this  introductory  stage 
the  "optimum"  material  possesses  the  greatest  potential  for  successful 
appl ications. 

Additional  fluidization  and  heat  transfer  studies  also  need  to  be 
conducted  using  varying  bed  and  heating  geometries.  The  results  of  this 
investigation  are  based  upon  tests  conducted  only  using  the  fluidized  beds 
and  heaters  described  in  Chapter  4.  These  results  have  then  been  assumed 
to  remain  constant  in  order  to  be  applied  to  the  solar  collector  and  heat 
exchanger  applications.  This  assumption,  although  reasonable,  must  be 
verified  through  experimental  means  in  order  to  determine  its  validity  and 
to  further  determine  optimum  conditions  for  low  density  particle 
fluidization  and  heat  transfer.  It  is  recommended,  for  example,  that  tests 
be  conducted  using  varying  finned  tube  geometries.  These  tests  would 
assist  in  understanding  the  effect  the  presence  of  fins  have  on 
fluidization  and  would  produce  further  data  available  for  selecting  an 
optimum  finned  tube  geometry  for  a  particular  application.  Similarly, 
additional  testing  on  the  effect  of  air  flow  obstructions,  roughened 
fluidized  bed  wall  surfaces,  and  varying  bed  depths  should  be  conducted  to 


aid  in  the  understanding  and  optimization  process  of  low  density  particle 
fluidization. 

An  additional  area  requiring  more  research,  which  is  of  particular 
importance  to  the  applications  discussed  in  this  work,  is  the  fluidized  bed 
thermal  resistance  determination.  The  theoretical  model  used  in  Chapter  7 
(Eq.  7.7)  is  based  upon  an  effective  heat  conduction  path  length  that  is  a 
function  of  immersed  heat  transfer  tube  spacing  (Eq.  7.8)  and  an  effective 
thermal  conductivity  (Eq.  7.9).  Recently,  though,  it  has  been  suggested 
that  a  more  realistic  model  expressing  this  resistance  should  be  based  upon 
a  path  length  that  is  a  function  of  the  average  bubble  diameter  rising 
through  the  bed  [41],  If  this  is  so,  it  can  be  expected  that  the  average 
bubble  length  may  be  less  than  the  length  represented  by  Eq.  (7.8)  and 
result  in  an  increased  bed  resistance  than  that  which  was  originally 
predicted.  The  magnitude  of  this  increase  will  determine  what  procedures 
should  be  taken  to  decrease  this  resistance  or  totally  re-evaluate  the 
feasibility  of  the  applications  considered  in  this  work. 

Consequently,  the  initial  successes  observed  in  this  investigation 
must  be  tempered  by  the  knowledge  that  many  more  questions  exist.  Not  until 
these  topics,  and  any  further  ones  they  in  turn  may  generate,  are  resolved 
will  a  full  understanding  of  low  density  particle  fluidization  and  heat 
transfer  with  successful  applications  be  realized. 
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APPENDIX  A 

WALL  HEATER  HEAT  BALANCE  AND  INSULATION  DIMENSIONS 


This  appendix  details  the  analysis  taken  in  determining  the  necessary 
insulation  required  to  be  placed  around  the  wall  heater  to  insure  that  the 
heat  flux  generated  is  transferred  into  the  two-dimensional  fluidized  bed. 
Figure  A-l  depicts  a  schematic  cross-sectional  view  of  the  wall  heater  in 
which  a  heat  balance  analysis  is  conducted. 

From  symmetry  and  the  geometry  of  the  heater  the  following 
equalities  exist: 


qtop 
qside  1 
Abed 
Atop 
Aside  1 


^bottom 
qside  2 

Aback  =  in*  x  ®  in*  =  *0^3  ft2 

Abottom  =  ^-5  in.  x  0.25  in.  =  .003  ft2 

Aside  2  =  ®  in*  x  ^*25  in*  =  ^2 


An  energy  balance  around  the  heater  yields 


qin  qbed  +  qback  +  qtop  +  ^bottom  +  2qside 

where 

qbed  ^bedAbed^wall  "  "'"bed^ 

a  =  k  A  (T  _  t  )/l 

^back  styrofoam  backv  wall  amb.'  styrofoam, back 

a  =  k  a  (T  -  T  )/l 

Mtop  styrofoam  topv  wall  amb.  styrofoam, top 


(A-l) 


[I] 


bed  bedv  wall 
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[  i  i  i )  _bed 


^bottom 


_ hbed  Sed^wall  ~  Tbedj _ 

k rubber  ^bottom ^wal  1  "  Smb.  ^bottom 


bottom 


TQQ(krubber)  Sottom^wall  Smb.^ 
hbed  Abed(Twall  "  Tbed} 


( i  v )  _bed 


hbed  Sed^wall  "  Sed^ 


^side  kstyr.  Side^wall  ”  "Smb.^Side 


=  100 


Side 


=  10°(kstvr.>  Side^wall  ~  Smb.^ 
Sed  Sed^wall  "  Sed^ 


Assumed  values  of  (Twall  -  Tamb< )/(Twa1 ,  -  Tbed)  =  5  , 
hbed  ■  20  Btu/hr-ft2-«F  ,  kstyrofoam  =  0.02  Btu/hr-ft-T  and 
krubber  =  Btu/hr-ft-°F  used  in  conjunction  with  Eqs.  (A- 
(A-5)  yield  the  required  insulation  thicknesses: 


Sack 

=  0.5  ft 

=  6.0  in 

top 

=  0.024  ft 

=  0.29  in 

l 

bottom 

=  0.11  ft 

=  1.29  in 

100 


(A— 4) 


(A— 5 ) 


through 


0.16  ft 


0.95  in 


APPENDIX  B 


FIN  ANALYSIS  FOR  FINNED  TUBE  HEATER 


This  appendix  details  the  particular  characteristics  of  the  fins 
used  on  the  finned  tube  heater.  These  fins  are  circumferential  fins  of 
rectangular  cross-section  and  the  procedure  used  in  analyzing  these  fins 
is  outlined  in  reference  [40]. 

As  illustrated  in  Fig.  4.5  each  fin  has  an  outer  radius,  xg  ,  of 
0.094  ft.,  an  inner  radius,  xfa  ,  of  0.031  ft,  and  a  thickness,  2yb  , 
of  0.0017  ft.  The  total  fin  area  available  for  heat  transfer  may  be  cal 
culated  using 

Afin  =  (No.  of  fins) (2  sides)  tt[x2  -  x2]  (B- 

=  13  ( 2 ) it  [(.094)2  -  ( .031  )2] 

=  0.64  ft2 


In  order  to  determine  the  fin  efficiency  the  term  lA\[Yy^  has  to 
be  determined  where: 


Z 

h 

k 
*b 


-  x. 
e  b 

10  Btu/hr-ft2-°F  (Assumed,  based  on  results  of 

unfinned  tube  data.) 

220  Btu/hr-ft-°F  for  copper  fins 

0.00085  ft. 


Using  these  values. 


2,/h/ky^  =  0.46 
and 

^e  =  .094  s 

.031 


(B-2) 

(B-3) 


yields,  from  reference  [40],  a  fin  efficiency  of  approximately  0.90. 

Further,  to  determine  the  finned  tube  heat  transfer  coefficient  the 
results  of  a  heat  balance  around  the  finned  tube  may  be  utilized 


qtotal 


qtotal 


%are  tube  +  qfin 
h  Abare  tube  ^  +  h  nfin  Afin 
h  AT[Abare  tube  +  (nA)fin] 


AT 


(B-4) 


assuming  that  the  heat  transfer  coefficient  of  the  bare  tube  and  the  fins 
is  approximately  equal.  From  Eq.  (B-4)  and  substitution  of  the  known  area 
and  calculated  efficiency  values  an  expression  for  the  heat  transfer 
coefficient  can  be  obtained 


qtotal 

(0.644)AT 


( B— 5 ) 


Equation  (B-5)  is  used  in  this  experiment  to  determine  the  finned 


tube  heat  transfer  coefficients. 


APPENDIX  C 


THE  THEORY  OF  STATIC  CHARGING 


In  Chapters  5  and  6  the  observance  of  electrostatic  charging  of 
the  Perlite  and  styrofoam  materials  in  contact  with  plexiglas  was  dis¬ 
cussed.  This  appendix  will  outline,  in  general  terms,  the  phenomena  of 
electrostatic  charging  between  two  materials  and  what  measures  can  be 
taken  to  reduce  this  effect.  This  discussion  will  aid  in  developing  a 
better  understanding  of  the  phenomena  as  observed  in  this  investigation 
and  for  future  work  with  similar  materials  in  fluidized  beds. 

The  subject  of  electrostatic  charging,  frictional  electricity,  or 
"triboelectric  charging"  has  a  history  of  having  been  long  recognized  but 
not  well  understood.  The  production  of  electric  charges  that  occurs 
when  two  solids  are  rubbed  together  (static  cling  of  clothing  for 
instance),  the  electrification  of  liquids  such  as  gasoline  when  they 
are  poured  out  of  a  container,  the  formation  of  charged  droplets  in 
spraying,  and  the  attraction  of  the  Perlite  and  styrofoam  materials  to  the 
plexiglas  bed  wall  are  all  examples  of  triboelectric  charging  of  materials 

It  has  long  been  recognized  that  if  two  dissimilar  metals  are  placed 
in  contact  there  will  be  a  potential  difference  between  the  metals.  This 
is  because  electrons  will  find  it  easier  to  go  from  one  metal  to  another. 
This  potential  difference  between  the  metals  is  called  the  contact  poten¬ 
tial.  If  the  two  metals  are  then  carefully  separated  one  will  be  positive 
ly  charged  and  the  other  negatively.  The  potential  difference  between 
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the  metals  will  be  at  most  3  to  4  volts,  but  is  usually  around  2  volts 
[35].  Similarly,  if  two  electrically  nonconducting  materials,  insulators 
or  dielectrics,  are  brought  together  they  may  also  become  charged.  If  a 
tangential  mechanical  force  acts  on  the  boundary  of  the  contacting 
materials,  the  force  must  overcome,  among  other  things,  the  electric 
attraction  between  the  two  dissimilar  materials  for  motion  occur.  As  a 
result  of  this  motion  the  two  materials  become  tangentially  separated.  A 
small,  but  finite,  fraction  of  the  total  mechanical  work  expended  is  spent 
in  producing  the  electrostatic  energy  of  the  separated  charges.  This 
electrostatic  energy  is  called  frictional  electricity.  It  has  generally 
been  accepted  that  the  process  of  material  charging  through  contact  and 
friction  with  another  material  is  all  considered  triboelectric  charging. 

In  general,  conductors  and  insulators,  when  electrically  isolated  from 
the  earth,  can  easily  become  charged  as  a  result  of  interacting  with 
another  material.  Friction  may  be  involved,  or  just  contact  followed  by 
separa  tion. 

Metals  can  transfer  electrons  from  one  surface  to  another  based  upon 
the  metal's  Fermi  energy  level.  When  two  dissimilar  metals  are  placed  in 
contact,  there  will  be  a  transfer  of  electrons  from  one  to  the  other  until 
the  potentials  of  the  metals  are  such  that  the  Fermi  levels  are  aligned. 

It  is  necessary  to  give  electrons  an  amount  of  energy  equal  to  the 
material's  work  function  in  order  for  them  to  escape  from  the  material's 
surface.  Consequently,  the  potential  difference  established  between  two 
metals  is  the  difference  in  work  functions.  Electrons  are  transferred 
through  the  interface  from  the  material  with  the  lower  work  function  to 
that  with  the  higher  work  function.  This  mechanism  of  charge 
redistribution  generally  holds  for  metal-to-metal ,  metal -to- semi  conductor 
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and  semi  conductor- to-semi conductor  contact  situations.  However,  for 
dielectrics,  such  as  Perlite,  styrofoam  and  plexiglas,  there  appears  to  be 
some  uncertainty  and  conflict  of  opinion  concerning  the  mechanism  of  their 
triboelectric  charging. 

Postnikov  has  postulated  that  in  contact  of  a  metal  and  a  dielectric 
charging  is  due  to  the  transfer  of  electrons  from  the  metal  to  the 
dielectric  and  the  transfer  of  positive  or  negative  ions  from  the 
dielectric  to  the  metal  [36].  Whereas  in  the  contact  of  two  dielectrics 
charging  results  from  the  diffusion  of  "charge  carriers"  from  one  substance 
to  another. 

Harper,  on  the  other  hand,  has  suggested  a  more  detailed  account  of 
the  static  charging  mechanism  of  dielectrics  [37,38].  He  agrees  that  the 
transfer  of  ions  plays  an  important  role  in  the  charging  of  solids,  but 
emphasizes  those  mobile  ions  that  have  been  formed  through  electrolytic 
dissociation  in  the  surface  films  of  moisture  on  solids.  He  suggests  that 
there  is  a  rupture  of  a  double  diffusional  layer  of  electrolytic  ions  in 
the  thin  water  film  between  two  contacting  dielectrics  or  between  a  metal 
and  a  dielectric  that  is  responsible  for  the  charging  of  the  surfaces  once 
they  are  parted  from  one  another. 

He  additionally  categorizes  dielectrics  into  two  classes.  The  first 
class  consists  of  "electrophilic"  dielectrics  which  acquire  strong  charges 


through  contact  with  metal.  These  include  glass,  fused  quartz,  and 
magnesium  oxide.  The  second  class  consists  of  "electrophobic"  dielectrics. 
These  materials  take  on  a  negligible  charge  when  in  contact  with  metal. 

This  class  includes  polyethylene,  polystyrene,  nylon  and  Perlite.  The 
observed  discharging  of  electrophilic  dielectrics  and  charging  of 
electrophobic  dielectrics  when  in  contact  with  metal  is  a  result  of 
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electrolytic  and  greasy  surface  contamination.  Harper  states  that  as 
greasy  (electrophobic)  contamination  is  accumulated  on  an  electrophilic 
dielectric  the  charge  will  decrease.  Similarly,  accumulation  of 
electrolytic  contamination  on  an  electrophobic  dielectric  will  increase  its 
charge.  He  further  considers  the  electrophilic/electrophobic  character  as 
an  intrinsic,  inherent  property  of  a  dielectric. 

This  conclusion  seems  well  founded.  According  to  classical 
electrostatic  theory,  the  greater  a  material's  ability  to  polarize  the 
greater  its  ability  to  be  charged  and  transfer  that  charge  [39].  The 
extent  to  which  a  dielectric  becomes  polarized  when  in  an  electric  field  is 
a  function  of  an  intrinsic  material  property  called  its  susceptibility,  X 
A  dielectric  coefficient,  K  ,  is  defined  as 

K  -  1  +  x  (C.l) 

The  larger  a  material's  susceptibility,  the  greater  is  its  ability  to  be 
polarized.  Consequently,  materials  with  a  larger  value  of  the  dielectric 
coefficient  will  be  polarized  easier  than  those  with  a  smaller  value. 

Table  C.l  is  a  listing  of  various  materials  and  their  corresponding 
dielectric  coefficient  value.  Note  that  water,  a  good  conductor,  has  a 
high  K  value  whereas  air,  a  poor  conductor,  has  a  low  one.  As  identified 
through  experimentation  by  Harper  and  shown  in  Table  C.l,  glass 
(electrophilic,  K  =  5-10)  is  more  capable  of  being  charged  when  brought  in 
contact  with  a  metal  than  polyethylene  (electrophobic,  K  =  2.25).  A 
similar  electrophobic  result  was  observed  in  this  investigation  with  the 
Perlite  and  styrofoam  contacting  the  metal  walls  and  heater  surface. 


TABLE  C.l  DIELECTRIC  COEFFICIENT,  K  [39] 


MATERIAL 

K 

VACUUM 

*1 .00 

AIR 

1.0006 

TEFLON 

2.1 

POLYETHYLENE 

2.25 

BENZENE 

2.28 

HEVEA  RUBBER 

2.94 

VINYLITE 

3.18 

PLEXIGLAS 

3.40 

MICA 

3-6 

BAKELITE 


5.50 
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At  this  point  it  would  appear  that  triboelectric  charging  results  from 
various  configurations  of  material  contact  charging  alone.  However,  it  has 
repeatedly  been  observed,  such  as  in  this  investigation,  that  when  certain 
materials  are  rubbed  together  they  become  electrically  charged.  Therefore, 
some  mechanism  must  account  for  this  and  two  hypotheses  are  commonly  used 
to  explai n  i  t. 

One  answer  to  the  rubbing  or  "friction"  question  that  has  been 
proposed  is  the  Vol ta-Helmhol tz  hypothesis  [37].  This  hypothesis  also 
attributes  all  triboelectrification  to  contact.  It  states  that  the  virtue 
of  rubbing,  however,  is  only  to  multiply  the  points  of  material  contact. 
When  two  materials  are  brought  together,  initial  contact  is  made.  When  the 
materials  rub  together,  old  contacts  are  broken  and  new  ones  are  made. 

Another  theory  accounting  for  frictional  electrification  is  through 
material  transfer  [37].  When  material  is  transferred  from  one  surface  to 
another  via  friction  the  material  will  be  removed  from  asperities  on  the 
surface  from  which  it  comes.  This  material  may  not  necessarily  be  typical 
of  the  surface  as  a  whole  and,  in  the  case  of  a  metal,  may  then  have  a 
different  contact  potential  from  the  rest  of  the  surface.  This  transferred 
material  could,  therefore,  carry  with  it  an  electric  charge  to  the  other 
surface.  If  the  material  transferred  comes  from  a  dielectric  hot  spot 
created  during  frictional  rubbing  i t  may  have  lost  electrons  to  the  less 
hot  material  left  behind  and  thus  carry  an  electric  charge  to  the  new 
surface . 

Thus,  the  mechanisms  of  triboelectrification  presently  understood 
entail  some  combination  of  contact  and  frictional  effects.  Depending  upon 
the  materials  involved,  their  geometry,  physical  and  electrical  properties, 
state  of  cleanliness,  nature  of  interaction,  etc.  one  or  both  of  these 
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processes  may  be  controlling.  Due  to  this  variety  of  variables  there  is  no 
universal,  single  criteria  or  relationship  for  determining  "a  priori"  how 
much  actual  electricity  is  generated  triboelectrically  under  known 
conditions  between  two  materials  in  any  situation.  Numerous  experiments 
have  been  conducted  in  search  of  this  relationship,  success  being  only  in 
developing  one  that  satisfies  the  particular  conditions,  materials,  and 
controls  of  that  experiment. 

However,  enough  is  known  about  this  phenomenon  in  order  to  minimize 
its  effect,  if  desired.  As  discussed  earlier,  electrically  conducting 
materials  have  a  large  dielectric  coefficient  value.  Thus,  metallic 
fluidized  particles,  for  example,  would  be  expected  to  exhibit  very  little 
static  electric  problems  in  fluidized  bed  operations.  Problems  would  be 
expected,  on  the  other  hand,  if  the  bed  walls  were  made  of  an 
electrophyll i c  material,  such  as  glass.  Similarly,  electrophobic  fluidized 
particles  would  be  expected  to  exhibit  little  static  electric  problems  in  a 
bed  made  with  metallic  walls.  In  this  investigation,  problems  occurred 
because  electrophobic  fluidized  particles.  Perlite  and  styrofoam,  were  in  a 
bed  that  had  an  electrophobic  material  wall,  plexiglas.  Triboelectric 
charges  were  built  up  on  the  dielectrics,  due  to  contact  and  friction,  that 
were  unable  to  be  sufficiently  transferred  away.  It  is  further  suspected 
that  the  styrofoam  has  a  lower  K  value  than  the  Perlite  since  it 
exhibited  the  strongest  triboe lectric  effect  with  the  plexiglas. 

In  some  engineering  or  laboratory  applications  it  may  be  impractical 
to  construct  a  fluidized  bed  or  specify  particular  particles  in  order  to 
reduce  the  triboelectric  attraction  phenomena.  In  these  instances, 
applying  the  theory  can  again  be  helpful.  For  example,  increasing  the 
relative  humidity  of  or  using  static  neutralizers  with  radioisotopes  in  the 


fluidizing  air  would  substantially  increase  its  dielectric  coefficient 
value.  Thus  the  fluidizing  air  may  become  conductive  enough  to  control  the 
accumulation  of  triboelectric  charges.  Adding  special  electrically 
conductive  fillers,  using  commercially  available  antistatic  preparations, 
or  simply  providing  an  electrical  path  to  drain  off  charges  on  the 
dielectrics  would  aid  in  increasing  their  electrical  conductivities. 

In  summary,  this  appendix  has  attempted  to  briefly  discuss  the  general 
aspects  of  triboelectrification  of  solid  materials.  The  mechanisms 
involved  in  this  process  are  a  combination  of  electrostatic  charging  due  to 
contact  and  friction  between  two  materials.  There  appears  to  be  no 
definite  agreement  on  the  exact  mechanism  that  explains  triboelectric 
charging  other  than  it  is  a  function  of  the  materials'  characteristics 
involved,  the  environment  in  which  they  are  contacted,  their  geometry  and 
the  nature  of  their  contact. 

Yet,  despite  this  lack  of  knowledge  of  the  exact  mechanism  involved, 
enough  is  known  to  satisfactorily  prevent  the  effect  of  triboelectric 
charging.  Additionally,  as  a  result  of  experimentation  and  experience  with 
triboelectricity  certain  general  relationships  or  "laws"  have  been 
developed  and  are  summarized  in  Table  C.2. 
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TABLE  C.2  SOME  "LAWS"  OF  TRIBOELECTRICITY  [36] 


SIGN  OF  CHARGE 

CONTACTING  MATERIALS  + _ ; _ 

Dielectric  and  fine  particles  Dielectric  Particles 

of  the  same  dielectric 


Dielectric  and  Dielectric 


<# 


Dielectric  with 
with  higher 
susceptibili ty 

Harder  Dielectric 


Lower 

susceptibil i ty 
Softer  Dielectric 


Metal  and  Dielectric 


Metal  and  Metal 


Dielectric 


Lower  Work 
Function 


Metal 


Higher  Work 
Function 


Softer  Metal  Harder  Metal 


l# 


APPENDIX  D 


SQUARE  EDGED  ORIFICE  DATA 

This  appendix  outlines  the  procedure  used  in  determining  the  air 
velocity  through  the  fluidized  beds  of  this  experiments.  Two  ASME  square 
edged  orifices  with  flange  taps  are  utilized  in  making  this  determination. 
The  specific  orifice  plate  data  presented  here  are  based  upon  the  general 
guidelines  listed  in  reference  [19]. 

The  equation  for  the  flow  rate  of  a  gas  through  the  ASME  square- 
edged  orifice  meter  with  flange  taps  is 


where 

a)  =  mass  flow  rate  [lbm/sec] 

D2  =  orifice  diameter  [in.] 

K  =  flow  coefficient 
Y  =  expansion  factor 

P-|  =  static  pressure  before  orifice  [in.  absolute] 

T.|  =  temperature  before  orifice  [°R] 

G  =  specific  gravity  of  gas  (air  =  1.00) 
y  =  supercompressibility  factor 
AP  =  pressure  drop  across  orifice  [in.  H^O] 


0.6056 


to  be  measured 
to  be  measured 


0.6050 


to  be  measured 
to  be  measured 


to  be  measured 


to  be  measured 


TABLE  D-l  Orifice  Parameter  Value  for  Flow  Rate  Equation 


Knowing  that  the  air  mass  flow  rate  through  the  orifice  can  also  be 


expressed  as 


w  -  P  UQ  A 


(D-2) 


where 


=  density  of  the  air  at  T-|  (p  =  0.075  lbm/ftJ  for  the  operating 
conditions  of  this  experiment) 

=  air  velocity 

=  cross-sectional  area  of  a  fluidized  bed 


yields  a  relationship  for  calculating  the  air  velocity  through  the 
fluidized  bed 


0.1145  K  Y  P. 

%  ’  < - TTT - >  <  t}  g  *  ap  > 


(D-3) 


Using  the  values  listed  in  Table  D-l ,  Eq.  (D-3)  can  be  modified  for 
orifice  A  and  orifice  B  use  in  the  two-dimensional  or  three-dimensional 
fluidized  bed.  Table  D-2  shows  these  modifications. 


ORIFICE  A 


ORIFICE  B 


P,  1/2  P.  1 

Two-Dimensional  Bed  J  3.733  (  =-  'P)  =  2.86  (  -J-  AP) 


APPENDIX  E 


TEST  MATERIAL  MINIMUM  FLUIDIZATION  VELOCITY  RESULTS 


This  appendix  presents,  in  both  graphical  and  tabular  form,  the 
results  of  the  test  material's  minimum  fluidization  velocity  determination 
experiments.  These  experiments  were  conducted  according  to  the  procedure 
outlined  in  Chapter  5.  The  average  Um^  for  each  particle  is  defined  as 
the  arithmetic  mean  of  the  two-dimensional  and  three-dimensional  fluidized 
bed  U  ^  values.  This  average  value  is  used  in  all  of  the  analysis  in 
Chapter  6  and  in  the  applications  of  Chapter  7.  Figures  E.l  through  E.7 
illustrate  the  experimental  results  and  Tables  E.l  through  E.7  list  the 
tabulated  results.  Also,  note  from  Figures  4.1  and  4.2  that  all  bed 
pressure  readings  are  based  on  a  pressure  tap  height  of  8  in.  for  the 
two-dimensional  bed  and  4  in.  for  the  three-dimensional  bed. 
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TABLE  E.1  NORTON  650  Umf  DATA 


APg  [in.  H20]  U  [fps] 


2- D  Bed:  0.74  0.12 

1.08  0.17 
1.50  0.25 
1.59  0.50 
1.58  0.77 
1.58  0.83 

3- D  Bed:  0.56  0.10 

0.64  0.10 
0.80  0.17 
0.84  0.72 
0.84  0.93 
0.84  1.11 


TABLE  E.4  NORTON 


APg  [in.  HgO] 


2-0  Bed:  0.46 

0.74 

1.02 

1.06 

1.07 

1.06 


3-D  Bed:  0.30 

0.40 
0.50 

n  co 


L«Ii] 
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TABLE  E.6  PERLITE  250  U  f 

APg  [in.  HgO] 

U  [fp 

2-D  Bed: 

0.16 

.005 

0.20 

.007 

0.25 

.009 

0.25 

.011 

0.25 

.014 

0.25 

.015 

3-D  Bed: 

0.045 

.007 

0.051 

.010 

0.095 

.014 

0.095 

.024 

0.095 

.031 

\1 


8-85 


DTIC 


